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SUMMARY 


This is the final report of a study of the applications of optical processing 
techniques and holography* The report includes studies of holographic informa- 
tion reduction, diffraction efficiency of photographically recorded interference 
patterns, recent advances in contour holography, hologram interferometry of 
transparent media, propagation through turbulent media, and the log-normal ran- 
dom process. 

Some preliminary analysis and experimental results of a new technique for 
reducing the information content of holographic records is presented. With this 
technique, three-dimensional images having a large viewing angle can be gener- 
ated at the expense of poorer image resolution. 

Some techniques for improving the diffraction efficiency of photographically 
recorded interference patterns are described. The results obtained with different 
bleaches and different bleaching procedures are reported. 

A new two-frequency contouring technique is described. With this technique, 
constant range or constant depth contours can be formed on image surfaces which 
are much larger than those previously obtainable. 

A preliminary study of hologram interferometry of transparent media is re- 
ported. The use of holographic interferometry to measure the index of refraction 
of homogeneous media has been explored, and an investigation of its application 
to heat and mass transfer and related fields has been initiated. 

In certain propagation experiments at optical wavelengths, it is well known 
that the first-order distribution of the received wavefront is log-normal. The 
basic way the log-normal law occurs is reviewed, and then an approximate solu- 
tion to the wave equation is obtained. A thorough discussion of the log-normal 
process is presented, including a derivation of its likelihood functional. 
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INTRODUCTION 

This is the final report on a study of optical data processing techniques and applications of 
hologram techniques. The study began on 1 March 1967 and* ended on 29 February 1968. This 
report includes primarily those contributions of the final six months of the grant period. The 
previous work performed under this grant has been reported [1] . 

A new technique for reducing the space -bandwidth requirements of a holographic recorder 
was described in the semiannual report [1] . The merit of this technique is that the information 
content of a hologram can be reduced radically without sacrificing the viewing angle. An analy- 
sis which indicates the improvement provided by this technique is presented here. Recent ex- 
perimental results support the analytical results. 

An investigation was undertaken both to examine old bleaches and to discover new ones 
which would improve the diffraction efficiency of photographically recorded interference pat- 
terns. Only three-dimensional diffraction gratings were considered. The diffraction of such 
gratings depends primarily on the index- of- refraction change which results from the bleaching 
process. Although diffraction efficiencies approaching 100% are theoretically possible, the 
highest measured value was approximately 47%. By eliminating surface reflections, this value 
could be increased to 60%. 

The generation of range contours holographically appears to be a very promising applica- 
tion of hologram techniques. Here, a modification of the previously reported [1] two-frequency 
technique is presented. With this new version, range contours can be generated over consider- 
ably larger surfaces than were previously possible. It is also shown that contour surfaces de- 
noting equidistant planes can easily be obtained. Such surfaces are generally much easier to 
interpret than ellipsoidal ones. 

Holographic interferometry has been used to measure the index of refraction of a trans- 
parent, homogeneous medium, and experimental results are illustrated. Also, its applications 
to heat and mass transfer have been considered. To illustrate the use of holographic interfer- 
ometry of nonhomogene ous media, we have exammed the density field of air surrounding heated 
objects. 

I The processing of holographically recorded images after they have propagated through 
turbulent media was first considered in the semiannual report. Here, the log-normal law rep- 
resentation of the received wavefront is reviewed. Then an approximate solution to the wave 
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equation is derived, definitions and properties pertaining to the log-normal process are pre- 
sented, and the likelihood functional is obtained. 
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HOLOGRAPHIC INFORMATION REDUCTION 


2.1. METHODS 

The need for reducing the information content of a hologram while maintaining, insofar as 
possible, some of the image characteristics (such as resolution,, viewing angle, field of view, 
or signal-to-noise ratio) has been well documented [2, 3, 4]. Three basic approaches to the 
problem have come to light. 

The first and most straightforward method is to perform a spatial filtering operation on 
the hologram, removing all but a narrow band of frequencies. For example, in the case of a 
hologram made and reilluminated with a collimated reference wave, the filtering operation re- 
duces the angular field of view which any small region of the hologram reconstructs. Thus, 
from any fixed position, an observer can view only part of the reconstructed image at a time, 
and he must scan the hologram to see the entire image. Such a system seems to have little 
use, although it might be of value in some particular application. 

Another method is to sample the complex object field in some plane between the object 
and the hologram. The number of samples recorded might vary from one or two samples [5] 
up to a very large number [6] , usually uniformly distributed. This method can reduce the 
spatial frequencies which must be recorded as well as the spatial bandwidth. However , the 
viewing angle is reduced, image resolution is slightly reduced, and usually a two-step process 
is necessary. In addition, the perspective becomes discontinuous, i.e., the image changes in 
discrete jumps as the observing position changes. The sampling method is useful when per- 
spective can be sacrificed rather than image resolution in order to gain a reduction in the in- 
formation which must be transmitted. 

The dispersion method [3, 4, 7] with which we are primarily concerned is fundamentally 
different. Full normal perspective is preserved while the reduction in information manifests 
itself in a reduced resolution in the reconstructed' image or in a decrease in the signal-to- 
noise ratio. A diagram illustrating this technique appears in figure 1. Light from the object 
is diffracted by a dispersion plane m(y) and falls upon the hologram plane h(z) where it is re- 
corded. If the hologram is constructed with a reference beam diverging from a point source 
and is then reconstructed with a beam originating from the opposite side of the hologram but 
which now converges to the same point, a conjugate wavefront traveling from right to left, 
again diffracted by the m(y) plane is generated. If the dispersion function m(y) is strictly a ' 
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FIGURE 1. CONSTRUCTION OF A HOLOGRAM THROUGH A DISPERSION PLANE 
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phase function, such as say a ground glass, then a pseudoscopic image may be observed to the 
left of the m(y) plane [8] . An observer viewing this image cannot see beyond Pj. The scatter- 
ing function appears to him to be the hologram, i.e., the source of the reconstructed wavefront. 
Thus, with a small hologram in Pg and a scattering function in P^, we have simulated a large 
hologram which appears to be in P^. The most important characteristic of this method is that 
it reconstructs an image over a large viewing angle from a relatively small hologram. 

2.2. THE INFORMATION-REDUCTION FACTOR 

To determine how much a given system reduces the information which must be recorded 
or transmitted, it is necessary to derive a standard means of comparison. The total informa- 
tion recorded by a one-dimensional plane hologram is proportional to the space-bandwidth 
product 

I = kLB (1) 

where k is a constant, L is the length of the hologram, and B is the spatial bandwidth {assumed 
here to be essentially constant over the entire hologram). For a hologram in plane P^ of fig- 
ure 2, with a point reference source S', and a field of view a, 

B = 2 sin a/2 (2) 

X 

Thus for L = 2Y * the information content is 
m 

k4Y sin a/2 

(3) 

Now, instead of recording the fringe system present in P^ on a photographic plate, let us re- 
place the plate, as shown in figure 3, with a scattering function m(y) whose extent is 2Y m , the 
same as the size of the previous hologram. Thus, the information is not recorded, but is 
still present within the aperture |y I ^ Y^. By a similar argument as above, the information 
recorded by a hologram in with reference source S 1 is 

j _ k2a sin p/2 ^ 

2 X 

The magnitude of 1 2 can be made arbitrarily small by decreasing the size of the hologram a or 
by increasing Increasing also reduces the maximum spatial frequency which the holo- 
gram must record, thus permitting the use of a lower resolution film or other detector. 

The reduction in information between a reference hologram recorded in and the actual 
hologram recorded in P 2 is 
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p, 

FIGURE 2. GEOMETRY FOR CALCULATING INFORMATION CONTENT OF A GENERAL 

HOLOGRAM' 
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FIGURE 3. GEOMETRY FOR CALCULATING INFORMATION CONTENT AT THE HOLOGRAM 

PLANE 
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(5) 


I, 2Y sin a/ 2 

rj _ _1 _ rn 

I ” a sin p/2 


In most applications, ~ Y^, « D^, and a « Y^. This allows the approximation 

sin P/2 

u 2 

which gives 


2D 0 sin a/2 

R = — - 

a 


( 6 ) 


In the case where X « D„ , we have 
m 1 


X 

sin a/2 s-=^ 

1 


thus permitting the approximation 


R s 


2X /D i 
m 1 

a/D n 


a 

e 


(7) 


We observe that in cases where these approximations are valid, the size of the scattering func- 
tion does not influence the information-reduction factor R. Since 6 can be arbitrarily smaller 
than very large reductions can be achieved. 


The above equations describe a one -dimensional model; for a square two-dimensional sys- 

2 

tern, the information reduction is R . 


2.3. ANALYTICAL DESCRIPTION OF THE SYSTEM 

It is impossible to achieve a reasonable degree of data reduction without sacrificing 
signal-to -noise ratio, resolution of the reconstructed image, or viewing angle. In order to 
examine what happens to the image, it is convenient first to present a general analytical de- 
scription of the system. 

The wavefront radiating from the object and falling upon the dispersion plane in figure 1 is 
g^(y), and that which is rediffracted by the dispersion plane and forms an image is g 2 (y). For 
simplicity in the analysis, the hologram aperture is assumed to be Gaussian weighted with a 
variable standard deviation of a. Using the conventional simplification of the FresneL diffraction 
formula, we may write the wavefront at the hologram plane as 
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where 

s = m(y)exp 

and F is the Fourier transform of s as given by 
s 

F g (u) = J s(x) exp (-j27iux) dx 

The total reconstructed image can then be considered to be a summation of such terms, 
given by equation 10, in which each term is representative of a point in the object space. 

In the following discussion, the characteristics of the image resulting from the use of a 
ground-glass dispersion structure are derived and demonstrated. It is shown that such a struc- 
ture causes the signal -to -noise ratio of the- image to be reduced as the recorded information is 
reduced (the hologram aperture is narrowed), A different type of random structure is then de- 
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scribed which causes the resolution of the Image to deteriorate rather than the background noise 
to build up. 

2*3*1* GROUND- GLASS DISPERSION FUNCTION. Assuming that the ground glass Is char- 
acterized by a random phase function, which has no transmission amplitude variations across 
the surface, choose a model in which the spatial bandwidth of m(y) is extremely broad and in 
which the phase relationship between spectral components is random. With an axial point ob- 
ject, the characteristics of the reconstructed image may be examined using equation 10. The 
width of the correlation function in (10) depends mainly on the frequency spread of the function 
FgtxADj} and has little dependence on the size of the hologram aperture* That Is, the resolu- 
tion of the point image depends on the extent of the m(y) function. However, the Intensity of the 
correlation peak is related to ’’a," the size of the hologram aperture. Because of the random 
nature of F s , the correlation function for off-axis points is relatively constant. As a becomes 
smaller, this correlation becomes larger. This Implies that the background noise in the image 
function is uniform and that it increases with decreasing aperture. In the limit, as a goes to 
zero, the term in the square braces is a delta function, and 


l*e,, the image function fg{x) is entirely noise. 

An experimental demonstration of these effects for a ground-glass dispersion function is 
shown in the image plane photographs of figure 4, Figure 4a is a photograph of the Image using 


these figures, it appears that reducing the hologram aperture has very little effect on the image 
resolution. However, there is a marked Increase in the background noise as predicted* 

2*3*2, RESOLUTION- REDUCING DISPERSION FUNCTION. In many cages it is not de- 
sirable to trade background noise in the image field for data reduction but rather to select a 
dispersion medium which causes the image resolution to decrease as H increases but which 
contributes very little to background noise, 

A dispersion medium which has this effect is again purely a phase function, but one with 
different statistics. The dispersion function can be written 


where 0 { y) is a random function of y. However, there are two essential restrictions on the 
Statistics of this random function. First, the derivative of the phase f>(y) must have a standard 




m(y) = exp [jflfyj] 
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FIGURE 4, IMAGING THROUGH GROUND -GLASS DISPERSION FUNCTION WITH INFORMATION CONTENT 






deviation which is large enough so that light from all parts of the object field Is diverted by the 
entire dispersion plane Into the relatively narrow hologram aperture. That is, it is essential 
that 


X + Y 


Y 

m m 


m 



D 2 


where a is the standard deviation of the phase rate L and 2X and 2Y are the dimensions of 

m m 

the object and dispersion function, respectively. When this requirement 3s satisfied, the large 

viewing angle is maintained. Next, the second derivative of the phase must be so smalt that G 

can be considered constant within a resolution element of the reconstructed wavefront* The 
\ 0 2 

size of this resolution element is — — . The second restriction implies that the phase relation- 
ship between spectral components is not entirely random; thus the previous discussion does 
not apply. Using a dispersion function having these special characteristics yields an image 
wavefront which is the conjugate of the object wavefront, except for an overlaying mask. 
When this mask is considered as a diffracting object, Its spatial frequencies may be directly 
related to a point- spread function in the image plane. For example, a component whose spatial 
frequency is I gives rise to an image-point location error of 


The equation relating the autocorrelation of the mask function to the autocorrelation of 6 
is (from sec. 2. 5) 




[ 0 1 - 


( 11 ) 


a 2 - kZ * 2 * o a 

P - •'ft' + 

i9 aut °correlation function of the mask function M, and R(t) is the autocorrelation 
function of b. 

From this equation, it is evident that the width of the autocorrelation function of R,.(t) de- 
bt 

creases as the size of the hologram aperture a decreases. The spatial frequencies of the Im- 
age spread increase, thus contributing to a poorer resolution. Similarly, as the correlation 
function for S becomes broader (that is. as fl becomes less rapidly varying), the image-point 
spread function narrows. However, there are practical limits to the breadth of Ffr), since a 
point is reached at which the Image will appear excessively spotty when viewed by eye. 
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From equation 11, 11 Is clear that the noise which can be attributed to the higher frequencies 

of R „(tJ does not build up as a is decreased. These desired characteristics are demonstrated 
M 

bv the experimental results shown in figure 5. The object was a ring surrounding a cross, 

2 

For figure 5a. the space- spatial bandwidth reduction ratio R was 140; for figure 5b, it was 500; 
and for figure 5c, it was 3600, Observe that the resolution does actually decrease in the manner 
predicted while the background noise is altered very little. 

In the actual experiment, the dispersion medium was an aluminized random surface on 
glass which, instead of transmitting light, reflected it into the hologram aperture. This was a 
convenient method for obtaining the large variation of diffraction angles across the surface. 

2.4, RELATED ACTIVITIES 

Some work has been done on producing a short animated movie by properly positioning a 
sequence of holograms in plane P 2 - The experimental problems encountered in such an at- 
tempt are very difficult, however, because of the precision with which these holograms must 
be located. Alternatively, a single photographic plate can be fixed in position in P,,, and the 
individual holograms of the movie can then be recorded side by side in some regular array. 

This is the method we used. The necessary scanning is performed by an opaque disc from 
which a 15° sector has been removed (fig, 6). This disc is driven by a stepper motor and its 
associated electronics. The stepping speed can be varied to produce different frame rates. 

A holographic movie which uses a reflecting scattering function is analogous to a conven- 
tional movie In that a wavefront is projected onto a screen. The distinction Ls that this projec- 
tion is coherent, and the image can be either in front of or behind the screen, or in the same 
plane. The advantage of a holographic movie made in conjunction with a scattering plate over 
a sequence of holograms alone is that small holograms can be used without reducing the viewing 
angle of the real image, 

2.5, DERIVATION OF EQUATION 11. 

In this section f an equation is developed which relates the correlation function of the overp- 
laying mask to the statistics of the dispersion function m(y). If the second derivative of the 
phase of the dispersion function is so small that the first derivative may be considered constant 
within a resolution element of the reconstructed dispersion plane, then 

y“ 

Hy) * i l \ 0 (u)du 

y - y j y 

a a E, VPfr> for |y’,y|<^ (12) 

y - y a 

Using this approximation in equation 9 and assuming an on-axis point object, one gets 
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FIGURE 5* IMAGING THROUGH A RESOLUTION -REDUCING DISPERSION FUNCTION 






FIGURE 6. SCANNING DISC ASSEMBLY FOR PRODUCING ANIMATED HOLOGRAM 

MOVIE 
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g 2 (y) ~ exp (— -^M(e) 


'1/ 


(13) 


where 


M(0) = J exp 


fe(y' 2 -2y’y)-^y' 2 

2D“ 


exp (jy’^dy’ 


The first exponential term in this equation is simply the point- image reconstruction func- 
tion. The function M is the mask function, and it is this portion of the above equation whose 
spectral density is examined. The problem of analyzing this spectral density, or alternatively 
the corresponding autocorrelation function, can be solved by using the transform method of 
nonlinear systems [9] , where 6 is assumed to be the Gaussian input to the system, and M(0) is 
the output. 


With this method, the autocorrelation R^(r) of the nonlinear system output is 

/ .2 2 \ 

~ Jf(w 1 ) exp 




2 2 \ / 2 2 \ 

U f H ^2\ 

— 2 — / <Ja, lJ F ^2^ exp l — 2 — ) dw 2 exp 


where F^) is the Fourier transform of M. That is, 


jk [ 2 _ \ -k 2 a 2 2 

2d\ 1 - 2w i y ).-^r w i 


F(ci^) = J M(0) exp (- d0 = exp 
Inserting this function into (14) gives 

k m w * | 


do). 


Jexp ^ 2 (|y 2 ) - 


( l ) 


2(T 
2 2 


du) r 


where 


°2 d 


After integrating, 


V 




R( T ) 


exp 


-k 2 / t 2 /3 2 2ry 2y 2 

2d 2 \/3 4 _R 2 /3 2 + R /3 2 + R v 


(14) 


(15) 


(16) 
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If there is to be any appreciable data reduction, either T)^ » or a « Y max * With 

these inequalities, the exponential term may be neglected. With the exception of a weighting 
function of y across the mask function, 


r m< t > - 


Vjs 4 - mr 


(id 


In an approximate manner, we can relate the standard deviation of the function to 

the image-point spread. For example, if the width of R m (t) is T, then the image-point spread 
is approximately AX - AD^/T. 


2.6. CONCLUSIONS 

The resolution- reducing dispersion function provided an image resolution loss without a 
decrease in signal-to-noise ratio. Another advantage of this type of dispersion function is that 
much less critical realignment is necessary during the reconstruction process than for the 
ground-glass system. A lateral misalignment on the order of the width of the correlation peak 
of R(r) will provide at least a blurry image. 

Although it is desirable to use as a dispersion function a surface which varies slowly, one 
which varies too slowly will cause the image to appear spotty. The size of these spots is re- 
lated to the width of the correlation slope function R(t). A compromise solution for the slope 
function must be chosen. 

In the dispersion- plane systems discussed so far, the pseudoscopic real image was recon- 
structed. Although this simplified the experimental work, the system is not limited to pseudo- 
scopic images. A different though related dispersion function can be used to reconstruct non- 
pseudoscopic -images. An example of a system which reconstructs an orthoscopic real image 
is one for which the actual image of the dispersion plane is focused by a telescope of unit mag- 
nification into a space which is readily accessible. In this focused plane, the conjugate of the 
dispersion function is placed. For the mirror dispersion-function system used to obtain the 
results of figure 5, this conjugate mirror is simply an imprint of the original surface. 

It should be mentioned that the dispersion function need not be confined to a planar function, 
but may use functions with depth or axial variations. An example of a system using a disper- 
sion function in two planes is one in which the first plane is an array of small diameter lenses 
which focuses the object onto a ground-glass plate in the second plane. This entire assembly 
may then be considered as the function m(y). Using such a system yielded results in which the 
resolution deteriorated with decreasing aperture. With this specific dispersion function, there 
may be some difficulty in avoiding focusing on the lens structure itself when the image is viewed 
by eye. However, it doqs illustrate the use of a three-dimensional dispersion function. 
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3 

TECHNIQUES FOB IMPROVING THE DIFFRACTION EFFICIENCY OF 
PHOTOGRAPHICALLY RECORDED INTERFERENCE PATTERNS 

3.1* INTRODUCTION 

Large interference patterns can be produced easily with lasers. Such patterns may be 
used as diffraction gratings, Fresnel zone plates, and other optical elements. Some uses al- 
ready have been described [10] , and other uses also exist. The major deficiency of such opti- 
cal elements is low efficiency, typically about 4%. With such an efficiency, possibly one ele- 
ment might be used in an optical system, but several in series would be impractical. 

An investigation was undertaken to find and examine bleaches which would improve the 
diffraction efficiency of photographically recorded interference patterns. If these patterns 
are recorded in thick recording media and if the density changes are converted into corres- 
ponding refraction index changes, then, theoretically, efficiencies of 100% could be achieved 
[11] . In practice, efficiency is limited by such factors as absorption m the emulsion, scatter- 
ing, surface reflections, and limited control over the index of refraction. With an ideal bleach, 
efficiency of 80% to 90% should be possible without the use of antireflection coatings on the 
surfaces. These efficiencies are for polarized light with the E- vector perpendicular to the 
plane of incidence. When the E- vector is parallel to the plane of incidence, diffraction efficiency 
is dependent on the angle of incidence. 

Experiments have shown that differenLbleaching techniques should be used for two-dimen- 
sional and three-dimensional diffraction gratings. For the two-dimensional case, bleaching 
techniques which produce relief patterns are efficient and easy to use, whereas for the three- 
dimensional case, the diffraction efficiency depends on index- of-refracti on change. Here we 
shall consider only the index- of- ref raction change. 

This investigation consisted of the following: (1) obtaining diffraction efficiency vs. density 
curves for various known bleaches, (2) finding the probable chemical reactions and identifying 
the chemical compounds causing the change in index of refraction, (3) investigating the losses 
in efficiency, (4) finding ways of stabilizing chemical compounds that tend to darken with ex- 
posure to light, and (5) finding other chemical compounds that might give the desired results. 

3.2. BLEACHING OF SILVER IMAGES 

To obtain the desired index- of- refraction change as a function of density, each exposed 
silver grain must be converted into a compound that has an index of refraction significantly 
different from that of the emulsion. The desired change can be expressed as follows: If we 
have a three-dimensional interference pattern described by T = k^ - kl (x, y, z), then we would 
like to replace it by a phase image such that after bleaching the transmission is one and the 
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index of refraction is n = n Q + k^x, y, z), where ^ is an appropriate constant. This gives 
phase variations in the emulsion of exp (yng) exp [yk^I(x, y, z)] . It is this type of phase vari- 
ation that we would like to achieve. 

In practice, each exposed silver grain is converted into a small dielectric particle having 
an index of refraction different than the emulsion. Variations in the concentration of these 
particles create an effective phase change. The maximum phase change should occur when 
the greatest difference exists between concentration of exposed silver grains in the light and 
dark fringe areas. If each silver atom would be replaced by one molecule of the bleach product, 
then all efficiency of diffraction versus density curves should have similar shapes. This is not 
the case, however, indicating that the mechanics of the bleaching process are more complicated. 

In the bleaching process, either the silver in the emulsion should react with the bleaching 
chemicals to form a compound that is insoluble in water or the bleaching solution, or the silver 
may act as a catalyst by which the bleach produces a compound at the location of the silver 
particles while the silver itself is eventually removed from the emulsion. With some bleaches, 
the silver both forms compounds in the emulsion and acts as a catalyst. The compounds formed 
in the emulsion should have the following properties: 

(1) The resulting compound should be transparent to light. 

(2) The refractive index of the compound should be significantly different from that of the 
gelatin. 

(3) The compound should be' stable and should not decompose or become opaque with ex- 
posure to light. 

All the tests for diffraction efficiencies were made using the recorded interference pattern 
between two spherical waves with an angular separation of 90°, as shown in figure 7. The spa- 
tial frequencies were about 2200 lines per millimeter. The interference patterns were recorded 
on Kodak 649F emulsion on 4 x 5 x 1/8-in. glass plates. The plates were developed in Kodak 
D-19 developer for 5 min, rinsed in short- stop bath for 15 to 20 sec, and fixed for 10 min. 
Continuous agitation was used in development and occasional agitation thereafter. The plates 
appear deep pink after fixing and become light pink after 20 min of washing in running water. 
After this processing, the plates were dried, the average transmission was measured, and 
then the plates were bleached. The detailed description of the bleaching process and of the 
probable chemical reactions is given in section 3.7. 

Figures 8 through 10 show the average transmission before bleaching versus efficiency- 
of- diffraction curves for various bleaches. The Kodak Chromium Intensifier Bleach (KCIB) 
gives the highest efficiency, 47% , and is the only one that is linear over most of its effective 
range. Various treatments were tried with plates processed in the KCIB to prevent them from 
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FIGURE 8. TRANSMISSION BEFORE BLEACHING VS. DIFFRACTION EFFICIENCY OF 

STANDARD BLEACHES 
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FIGURE 9. TRANSMISSION BEFORE BLEACHING VS. DIFFRACTION EFFICIENCY OF.R-IO 

SERIES BLEACHES 
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FIGURE 10. TRANSMISSION BEFORE BLEACHING VS. DIFFRACTION EFFICIENCY OF 

SPECIAL BLEACHES 
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turning dark after processing* If the plate is soaked in mercuric chloride (HgCl^) bleach for 5 
min as a final step and not washed afterwards, the plate remains clear* Because of the ex- 
tremely poisonous nature of HgClg, this is not a desirable treatment. We found that soaking in 
copper bromide (CuBr 2 ) or copper chloride (CuC^) also prevents the plate from turning dark. 

Besides the chromate-type bleaches, several others were found to improve the diffraction 
efficiency. Among these the potassium ferrocyanide (K^Fe(CN) 6 ), copper bromide and hydrogen 
peroxide (CuBr 2 + potassium bromide with potassium ferrocyanide (KBr + K^Fe(CN)g) 

were some that were quite good. The first two do not cause the plate to become more opaque 
with exposure to light or air, and the KgFe(CN)g by itself produces a very stable compound. 

The maximum measured diffraction efficiencies of the above three compounds were 39%, 37%, 
and 34%, respectively. The bleach containing hydrogen peroxide is difficult to use because 
oxygen bubbles tend to form inside the emulsion and tear it away from the base. This effect 
can be reduced if a low concentration of is used (1% solution) and bleaching time does 
not exceed 10 min, 

3.3. MEASUREMENT OF DIFFRACTION EFFICIENCY 

The diffraction efficiency as used here is defined to be I .. £t /l. , where I. is the intensity 

J diff me’ me J 

of the incident beam of light and I^. ff is the diffracted intensity of light. The plates were re- 
placed in the same position as when exposed (as shown in fig. 7), with one beam blocked and the 
other collimated. Collimation was used to eliminate dependence of the light meter on distance 
from the plate and to obtain a more intense beam of light. A YSI-Kettering model 65 radiometer 
was used to measure light intensity. The plates were always aligned for maximum diffraction 
efficiency and sometimes required a few degrees rotation relative to their original positions. 

The readings were corrected for changes in beam width due to rotation of the plate. 

In general, it was difficult to obtain repeatable results for a given exposure and process. 

The main problem seems to be the distortion of the emulsion caused by stresses in the emul- 
sion or the processing. The curves include losses due to surface reflections which could be 
eliminated by coating the surfaces with antireflection coatings. To obtain diffraction efficiencies 
with reflection losses removed, the diffraction efficiencies should be multiplied by a factor of 
1.23. 

‘ 3.4. LIGHT DISTRIBUTION IN A BLEACHED GRATING 

To make further improvements in bleached diffraction gratings, it is important to know the 
losses that occur in a grating. One grating, bleached in modified R-10 bleach and having a dif- 
fraction efficiency of 44%, was examined for losses. The losses were assumed to be caused 
by absorption in the emulsion, surface reflection, scattering, and transmission through the 
plate. 
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The scattered light was estimated by measuring the light intensity at a large distance — 

compared to the size of the grating — on either side of the plate and perpendicular to it, taking 

the average, and assuming the scattering to be uniform. If I T is the total scattered light inten- 

1 47tD^ 

sity and I is the average of two measured values, then I T - I M A - , where D is the distance 
from the grating and A is the area of the detector. The surface reflection from a single air- 
glass surface was measured at 9%, and at the second surface the energy loss is estimated at 4‘ 


The diffracted light and transmitted light were measured directly, and the difference of the 
measured intensities and the input is assumed to be absorbed. The light distribution for the 
sample measured is as follows : 


Diffraction 

44% 

Transmitted light 

8% 

Surface reflections 

13% 

Scattered light 

9% 

Absorbed light 

26% 


3.5. CHEMICAL REACTIONS IN BLEACHING 

The chemical compounds that are produced in the bleaching process are of interest, since 
the knowledge of these helps in the understanding of the bleaching process and helps in finding 
bleaches that have more desirable properties. Some chemical reactions were found in the lit- 
erature, while others were assumed to take place but may not be correct. Several of the 
bleaches are discussed here. 

The mercuric chloride bleaching process has been described by R. W. Green [12]. The 
reaction was described as Ag + HgCl 2 - AgHgCl^xHgO. The plate has a tendency to turn dark 
with exposure to light, indicating that it decomposes into free silver or mercury. If some 
HgCl 2 is left in the emulsion, the bleach is more stable. 

Bleaches containing chromates gave the best results, and these were studied in more de- 
tail. Pure silver particles were treated with various bleaches, but only the KCIB reacted with 
the silver in bulk form. The characteristic brownish color was not produced and the only pre- 
cipitated product seemed to be silver chloride (AgCl). The clearing bath (part B of KCIB) of 
the Kodak bleach did not react with this precipitate and apparently reacts with other deposits 
in the emulsion. 

One reaction produced by bleaching with a chromate bleach [13] is 

6 Ag + K 2 Cr 2 0 ? + 7H 2 S0 4 - 3 Ag 2 S0 4 + Cr^SO^I + K 2 S0 4 + 7H 2 0 

Only chromium sulphate (Cr 2 (S0 4 ) 3 ) remains in the emulsion. The diffraction efficiency of this 
bleach is low. This is probably because of the low refractive index of Cr 2 (S0 4 ) 3 , but it also 
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may be because of the low ratio of chromium molecules to silver atoms (1:6), and possibly the 
formation of Cr^SO^ molecules away from silver atoms. 

It appears that most of the good bleaches contain chlorine or bromide to form AgCl or 
silver bromide (AgBr) in the emulsion, which usually give high diffraction efficiencies. In the 
Kodak R-10 bleach, the probable reaction products remaining in the emulsion are Cr^SO^ 
and AgCl, but the AgCl seems to have more effect. The most efficient chromate bleaches ap- 
parently leave AgCl and chromium chloride (CrCl^) in the emulsion. 

Two bleaches were found that do not contain chromium or chloride and still have rather 
high diffraction efficiency. One of these is potassium ferrocyanide K 3 Fe(CN)g, which probably 
forms either AgK^Fe(CN)g or Ag^Fe(CN)^ in the bleaching process. The resulting compound 
is very stable, reasonably transparent at red wavelengths, and gives a maximum efficiency of 
39%. The other bleach contained copper bromide CuBr^ and hydrogen peroxide H 2 0 2 , which 
is believed to form two compounds, AgBr and copper oxide (Cu^O) in the emulsion. Both of 
these have high indexes of refraction, 2.25 and 2.70, respectively. The bleached plate has a 
light brownish color. It can be further processed by washing it in KCIB clearing bath, which 
turns it white, and then soaked in HgCl 2 to stabilize the compound. Since the composition of 
the clearing bath is unknown, the final compounds in the bleached plate are also unknown. Many 
other chemical compounds are known to have high indexes of refraction, some of which are 
listed in table I. 

'3.6. DISCUSSION 

The conversion of silver grains into transparent salt is a process dependent on many fac- 
tors not well understood. The efficiency of the bleached grating likewise depends on many fac- 
tors, including the handling of the exposed plates in the bleaching process, stability of optical 
systems during exposure, relative amplitude of the two interfering beams, and distortions in 
the emulsion. Because of these factors, it is difficult to repeat experiments, and slight vari- 
ations in efficiencies usually exist between gratings prepared in a similar manner. The graphs 
were drawn using the best efficiencies for each exposure. 

None of the bleaches that we have investigated possesses all the desired properties. The 
two most serious shortcomings are that the bleached plates still absorb light and that they tend 
to turn dark with exposure to light. Fixing apparently removes compounds responsible for 
index- of- refraction changes and, therefore, plates are not fixed after bleaching. The index- of- 
refr action changes also are not sufficiently high. Treatment with HgCl solution after bleaching 
and clearing seems to prevent darkening of the gratings. 

Other bleaches were tried. Good results were obtained with Kodak In- 4 intensifier bleach. 
Up to an optical density of 2.8, it follows the Kodak bleach, at which point it flattens out with a 
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TABLE I. CHEMICAL COMPOUNDS WITH A HIGH INDEX OF REFRACTION 


Compound 

Index of 
Refraction 

Color of Bulk 
Compound 

AgCl 

2.071 

White 

Agl 

2.21 

Yellow 

AgBr 

2.253 

Pale yellow 

HgO 

2.37 

Yellow or red 

Hg 2 CI 2 

1.97 

White 

Hg 2 Cl 2 • 2H 2 0 

2.35 


HgCl 2 

1.859 

Colorless 

CuO 

2.63 


Cu 2 0 

2.705 

Red 

Cu 2 0 3 

2.5 

Green 

CuBr 

2.160 

White 

Hgl 

2.75 

Red 

SrS 

2.11 

Colorless or light gray 

T1C1 

2.25 

White 

T1 I 

2.78 

Yellow or red 

Snl 4 

2.106 

Orange -red 

A1 4°3 

2.7 

Yellow-green 


2.78 

Ruby-red 

BaO 

1.98 

Colorless or white- 
yellowish powder 

BaS 

2.16 

Colorless 

CdO 

2.49 ’ 

Brown 

CdS 

2.506 

Yellow-orange 

Cr 2 (S0 4 ) 

Unknown 

Unknown 

CrCIg'xHgO 

Unknown 

Unknown 
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maximum efficiency of about 33%. Another bleach (not. standard) used KBr and K^Fe(CN)g and 
gave an efficiency curve that behaved like modified R-10 with KBr up to an efficiency of 30%, 
at which point it leveled out, reaching a maximum of 35%. 

Different developers were also tried: D-19 and BRP gave equally good results, while 
Dektal and D-8 showed a decline in efficiencies* The results given here are for a Kodak D-19 
developer. 

Several bleaches are available which would give diffraction efficiencies in the range of 35% 
to 48%. This variable represents about one-half f-stop change and, in ordinary photography, 
would be considered insignificant. Within this range, diffraction gratings or similar optical 
elements could be made quite easily. Some difficulty is expected in obtaining uniform diffrac- 
tion efficiencies over larger areas of film since the emulsion tends to get distorted. Without 
surface reflections, the best efficiency achieved would be 60%. 

The many problems and unanswered questions concerning bleaching provide a challenging 
project for a chemist or a thesis project for a chemistry student. 

3.7. CHEMICAL REACTIONS AND BLEACHING PROCEDURES 

3.7.1. KODAK CHROMIUM INTENSIFIER BLEACH 

Chemical Reaction: Ag + + 4C1 + Cr +++ -*■ AgCl \ + CrCl^ I 

The chemical compositions of the bleach and the clearing bath are unknown; therefore, 
the chemical reaction is not definitely known. 

Procedure : 

(1) Mix bleach and clearing bath according to manufacturer’s instructions. 

(2) Bleach for 1 min after the silver image has been removed (approximately 3 to 5 min). 
Agitate. 

(3) Rinse 0.5 min in running water. 

(4) Soak in clearing bath (part B) .until exposed emulsion becomes white, for about 5 to 20 
min. Use occasional agitation. 

(5) Wash in running water 5 to 10 min. 

(6) Soak 5 min in Photo- Flo and CuCl^ solution. 

(7) Dry without washing. 

3.7.2. KODAK R-10 BLEACH AND MODIFIED R-10 BLEACHES. 

Chemical Reactions: 
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(1) Standard R-10 bleach with NaCl in solution B: 

6Ag + 7H 2 S0 4 + 12NaCl + (NH 4 ) 2 Cr 2 0 ^ 

6AgCl \ + 2CrCl 3 * + 7H 2 0 + 6Na 2 S0 4 + (NH 4 ) 2 S0 4 

(2) Modified R-10 bleach with KBr in solution B: 

6Ag + 7H 2 S0 4 + 12KBr + (NH^Cr^- 

6AgBr \ + 2CrBr 3 * + 7H 2 0 + 6K 2 S0 4 + (NH 4 ) 2 S0 4 

(3) Modified R-10 bleach with KI in solution B: 

6Ag + 7H 2 S0 4 + 12KI + (NH 4 ) 2 CR 2 0 7 - 

6AgI| + 2CrI 3 * + 7H 2 0 + 6K 2 S0 4 + (NH 4 ) 2 S0 4 


Composition of bleaching solutions [14] : 

Stock Solution A 

Distilled water 500 ml 

Ammonium bichromate 20 g 

Concentrated sulfuric acid 14 ml 

Distilled water to make 1000 ml 

Stock Solution B 

Sodium chloride 45 g 

Distilled water to make 1000 ml 


For the bromide or iodate modifications, 92 g potassium bromide or 128 g potassium iodate 
instead of the sodium chloride was used in solution B. The NaCl in solution B was also replaced 
with NaF, but the effect was to remove emulsion without bleaching the silver. For use, mix 1 
part A to 1 part B to 10 parts water. 

Procedure: 

(1) Bleach for 1 min more than required for the removal of the silver (3 to 5 min). Agitate. 

(2) Rinse in running water, soak 1 min in CB-6 (clearing bath), and wash (approximately 
5 min). 

(3) Soak for 5 min in respective cupric halogen (CuClg, CuBr 2 , etc.), which also has Photo- 
Flo in it. The cupric halogen solutions contained 5 g of cupric halogen per 100 ml 

of water. 

(4) Dry without washing. 

3.7.3. POTASSIUM FERRXCYANIDE BLEACH 

Chemical Reaction: Agj + K 3 Fe(CN) 6 - AgK 3 Fe(CN) 6 1 or 

4Ag | + 4K 3 Fe(CN) 6 -Ag 4 Fe(CN) 6 * + 3K 4 Fe(CN) 0 
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Procedure: 

(1) Bleach for 4 min in 15 g of K 3 Fe(CN) 6 per 1000 ml of water. 

(2) Wash for 15 min in running water. 

(3) Soak in print flattening solution for 5 min. 

3.7.4. COPPER BROMIDE BLEACH 

Chemical Reaction: 4Ag j+ CuBr 2 + H 2 0 2 - 4AgBrj + CugO* + H 2 0 

Note: the copper may also form CuO or CuOg salt. 

Procedure: 

(1) Mix a solution in the proportions of 40 ml of 15 g of CuBr^ per 1000 ml of water solu- 
tion with 60 ml of diluted hydrogen peroxide (1 part of 3% hydrogen peroxide with 2 
parts of water). Mix the two solutions immediately before use because 0 2 is liberated 
from the mixed solution: 2H 2 0 2 — 2H 2 0 + 0 2 in presence of CuBr 2 . 

(2) Bleach the plate until exposed emulsion turns brownish, but not longer than 10 min. 
Continuously brush off bubbles forming on top of the emulsion. 

(3) Wash plate for 15 min in running water. 

(4) Rinse in Photo- Flow and dry emulsion, or continue with the following steps which 
make the plate more transparent and improve the efficiency somewhat. 

(5) Soak the plate for 8 min in KCIB, part B, the clearing bath. 

(6) Wash for 5 min. 

(7) Soak the plate in HgCl solution for 2 min. 

(8) Dry without washing. 

Note: The HgOg has a tendency to decompose thus forming oxygen bubbles inside the emulsion, 
and, in some cases, the emulsion will peal off. For this reason, the concentration of the bleach 
should not be stronger than indicated above, and the bleaching time should-not exceed 10 min. 

3.7.5. MERCURIC CHLORIDE BLEACH 
Chemical Reaction: 

2Ag + 2HgCI 2 - 2AgCH + Hg 2 ++ + 2Cl" 

Ag + HgClg - AgHgCl 2 xH 2 0 
Procedure: 

(1) Bleach 3 to 5 min or 1 min after dark silver image has disappeared. 

(2) Wash for 4 min, Photo- Flo rinse, and dry. 
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4 

RECENT DEVELOPMENTS IN CONTOUR HOLOGRAPHY 
4.1. INTRODUCTION 

The two-frequency method has been used successfully to generate constant range contours 
on or near holographic images [15]. One such method for obtaining the desired contours is to 
construct a hologram of an object at one wavelength then the hologram, as well as the ob- 
ject, is illuminated with a different wavelength A ^ By carefully positioning the holographic 
image in near coincidence with the object, an interference pattern results which is comprised 
of contours denoting constant depth. One problem with the above procedure is that the holo- 
graphic image becomes translated, rotated, and distorted when the reconstruction wavelength 
differs from the construction wavelength. Consequently, contours can often be produced only 
on a relatively small portion of the object surface. However, this problem can be greatly al- 
leviated by appropriately modifying the geometry for constructing the required hologram. 


Consider a small portion of a hologram constructed with a collimated reference of wave- 
length X 1 and then interrogated with a collimated beam of wavelength Each point of the 
object at a distance R 1 from a small hologram is imaged at a distance R^ such that 



(16) 


Thus, to minimize the relative change from the object to the corresponding image position, 
the object- to -hologram distance should be as small as possible. If plane- wave reference and 
readout beams are used, this condition can be relaxed considerably by utilizing only the object- 
to-hologram rays which are nearly longitudinal. Then the displacement of the image of the 
hologram is in a direction having the least effect on the range contours. 

We shall now present a practical method for generating range contours over large surfaces. 
Often one to two orders of magnitude improvement can easily be obtained for a given wavelength 
ratio. Furthermore, by properly using both collimated reference and illuminating beams, we 
can obtain contour surfaces denoting equidistant planes parallel to the hologram.* These sur- 
faces are generally much easier to interpret than the ellipsoids obtained by previous techniques. 


*Plane contour surfaces can be obtained using techniques previously reported, for ex- 
ample, by requiring all beams to be collimated and then viewing the reconstruction with a 
telescope. However, the maximum object slope over which contours can be thus recorded 
is considerably less than that obtainable with techniques described below. 
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4,2. THEORY 

Rather than attempting to place the object very near to a hologram (e.g., by illuminating 
the object from the back side of a hologram), we image an exact replica of the object near the 
hologram by utilizing an imaging system. Consider the unit telescope system shown in figure 
11, This telescope performs two basic functions. First, it produces an image I of the object 
O such that a hologram can easily be constructed in a plane which passes near or actually 

through the image. Second, its aperture A restricts the image-forming rays to be nearly lon- 
gitudinal. 


To obtain contour fringes of the diffuse object 0, we first construct a hologram H located 
very near to the image I with wavelength X^. Then the hologram and the object are illuminated 
with a new wavelength -By carefully positioning the image of the hologram in near coinci- 
dence with the image of the telescope, an interference pattern results denoting constant range 
contours. If a collimated reference beam is used to construct the hologram, the appropriate 
reconstructing beam is also a collimated'beam but of slightly different orientation than the 
reference beam. In particular, we require that the directions of the reconstructing beam a 
and the reference beam a satisfy the following relationship: 


X 

sin a = v 
c X 


2 

1 


sin a 

r 


(17) 


Let us now consider minimum requirements for the formation of contour fringes on or 
near the image of the hologram. It is well known that contour fringes occur due to interference 
effects between two images. Therefore, the lateral and longitudinal displacements of the re- 
constructed image should be less than the resolution and depth of focus of the viewing system, 
respectively. However, we restrict either thedmaging system or the subsequent viewing sys- 
tem to small viewing angles relative to their axes. Thus, no appreciable lateral displacements 
result, and only longitudinal displacements need to be considered. 

The resolution, depth of focus, and viewing angle are, respectively, 

X N f 

p “ 9 ~ *d ( 18 ) 


<»> 

_ d 

® = f (20) 

where \ denotes the wavelength, f the focal length of each lens, and d the diameter of the aper- 
ture A. According to equation 16, the longitudinal displacement of the image of the hologram is 
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( 21 ) 


where 


X 1 ' X 2 


AX = X^ - X 2 


X + X 


The above discussion on interference requires the following necessary condition: 


z i ■ z i 


After substituting (19) and (21) into (22), we find the following constraint on the longitudinal 
distance from the hologram to its reconstructed image 

l_ I c 4X 2 4X 2 (i \ 2 


0 2 Tax! \d, 


Of course, if the entire image is to be recorded, its longitudinal extent must be less than the 
depth of focus associated with the recorder. Thus, we usually need the more stringent con- 
dition 


M = ^ =2x/ ® 2 = 2 x ( h ) (24) 

Let us now consider the light intensity which results when the image of the hologram and 

that of the imaging system are made to interfere. Assume that the two image intensities are 

equal. It is convenient to define a rectangular coordinate system having its z-axis normal to 

the hologram and the elements of the imaging system while the x- and y-axes are in the plane 

of the hologram. A point source located at (x , O, z ) illuminates the object. The desired im- 

s s 

age intensity of an arbitrary object point (x 0 , y Q , z Q ) is then given by 

I' = 21 (l + .b cos jo^ - k 2 )z Q + (kj - k 2 ) 

[(z 0 - zf + (x 0 - x s ) 2 + y 2 ] 7 + /sjj (25) 

where I denotes the image intensity in the absence of the reconstruction of the hologram; b is 
some real constant between 0 and 1, depending upon the amount of image displacement associated 
with the holographic image; and are the wave numbers corresponding to X^ and X and /3 
is an unknown constant phase. According to (25), the contour maxima occur when 
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for integer values of n. 

Equation 26 denotes a set of paraboloids of revolution about the line (x - x g , y = 0) with 

vertices at ^x g , O, However, we are interested* in generating plane 

contour surfaces. Such surfaces can easily be obtained by merely letting \z I get very large. 

s 

More specifically, equation 26 becomes 

z n = noTrri + constant (27) 

0 2 |AX| v 


provided that 

(2|x s l + l x 0 1 ) 2 + Yp ^ V 2 

l z s - z q I ” 


where the lateral extent of the object is bounded by 


and 
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(28) 


Inequality 28 was obtained by requiring the resulting phase error to be bounded by ±tt / 4 
radians. Plane contour surfaces can be obtained with the setup shown in figure 11 by merely 
replacing the point-source illuminator with a collimated one oriented along the z-axis. This 
can easily be done .by incorporating a beam splitter between the object Oand the objective lens 
to reflect a collimated wave onto the object. 

As the angle between the normals of the hologram and the object is increased, the fre- 
quency of the contour fringes increases. The object orientation for which the contour fringes 
become indistinguishable will now be determined. Referring to the object shown in figure 12, 
we observe that separate contour fringes are imaged provided that 
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Combining (23) and (29) to eliminate 8 yields 
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This important result shows that contours can, in theory, be generated over surfaces of ar- 
bitrary large lateral extent. The maximum allowable inclination of the object is determined 
by the position of the object as well as the two wavelengths which are used. 

Practical limitations actually require that the right side of inequality 30 be greatly re- 
duced, As previously noted, if the entire object is to be recorded, the allowable range of z 1 
values must be reduced by a factor | AA |/A. Furthermore, if good fringe contrast is to be ob- 
tained, still greater reduction in the extent is required. Consequently, inequality 28 becomes 
after practical considerations 

|z_l f (tan = o (31) 

1 2 (A 

- Clearly, there is nothing unique about the imaging system shown in figure 11. In fact, a 
lensless system can be implemented by employing holographic techniques. Consider a doubly 
exposed hologram constructed by sequentially exposing a photographic plate to two different 
wavelengths, say A^ and A^. Rotating either the reference beam or the hologram between ex- 
posures allows us to read out the images separately. First, the exposure made with wave- 
length is read out with a A^ beam to form a real image in the original object space. This 
image is then recorded by a second hologram situated in the original object position. Next, 
the second exposure of the first hologram is read out with' wavelength Ag. Simultaneously, the 
hologram constructed in the object plane is read out with the appropriate A 2 beam to form the 
desired contour fringes. Clearly, such a doubly exposed hologram duplicates the task per- 
formed by the unit telescope of figure 11, except for the aperture A which restricts the image- 
forming rays to being nearly longitudinal. However, such an aperture can be incorporated into 
the subsequent recording or viewing system. It should be noted that the holographic imaging 
system just described causes two pseudoscopic rather than orthoscopic images to interfere, 
but no noticeable change in the recorded image and accompanying depth -contours results. 

The final contouring system which we consider requires one large lens and one large pho- 
tographic plate. Consider two small holograms constructed at the aperture plane A of fig- 
ure 11, with respective wavelengths A 1 and A 2 . Now by illuminating the two holograms with 
beams which are conjugate to the original constructing reference beams, images at both wave- 
lengths A^ and A 2 can be generated at the original object position. One of these images can be 
used to construct a third hologram in a plane near the original position. The image of this 
third hologram can then be made to interfere with the image of the other hologram recorded 
at plane A. In essence, the system just described uses one lens for both the. objective and the 
eyepiece lenses of a unit telescope, except for the fact that pseudoscopic images are formed. 
The size of the frequency-plane aperture A is determined by the extent of the holograms con- 
structed there. 
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FIGURE 13- DEPTH CONTOURS ON FOUR FLAT SURFACES 
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4*3, E XP E RIMENTAL RESU LTS 

The system shown in figure 11 has been used to generate constant depth contours* A 
Spectra Physics model 135 helium neon (He-Ne) laser capable of generating both a red line 
(X^ = 6 32 BA I and an orange line fx^ - GllflA) was used. Thus, consecutive contour fringes de- 
note a depth separation of 

* 1*2 

nnuT = 025 

The object used is comprised of four flat surfaces* each of slightly different angular orienta- 
tion, and a fifty -cent coin. The lateral extent of the object on which we were able to record 
fringes was about 10 cm by 0 cm. Actually, the lateral dimensions over which we could ob- 
tain depth contours— as the above theory predicts — was only limited by the useful apertures 
of the lenses used to image the object and those required to collimate the illuminating and 
reference beams* 

Figure 13 is a result obtained by using a telescope having a dc magnification of 2 and a 
viewing angle 0 of about 1/50 rad. The depth contours can be seen on the four flat surfaces 
as well as the coin. The curvature associated with the contours located at the extremities of 
the photograph is primarily attributed to aberrations introduced by our collimating lens. 

4.4. CONCLUSIONS 

We liave shown that two -frequency range (depth) contours can be formed over areas con- 
siderably larger than were previously possible. This is accomplished by using an imaging 
system so that the required hologram can be constructed arbitrarily close to the resulting 
Image. Additional improvement Is obtained by using only the imaging rays which are nearly 
longitudinal. Working with plane-wave sources, it is possible to form contour surfaces which 
are plane and equidistant from the hologram. 

The contouring techniques which we have described are not limited to optical wavelengths. 
Clearly other electromagnetic wavelengths, as well as acoustical sources, could be used. 
Generally, the required imaging system does not have to be of very high quality, providing 
that It functions the same way at both the wavelengths which are used. The basic result given 
In this report Is that the Lateral extent over which range contours can be obtained is only 
Limited by the longitudinal extent of the object surface and by its maximum slope. 


5 

HOLOGRAM INTERFEROMETRY OF TRANSPARENT MEDIA 
5.1* INTRODUCTION 

We have investigated the use of holographic interferometry to study phenomena related to 
the measurement of refractive Index of transparent media. The principle underlying this 
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important application of holography is that the wavefront emanating from a coherently illumi- 
nated object can be holographically reconstructed so precisely that it is optically indistinguish- 
able from the actual wavefront from the object. The actual and reconstructed wavefronts can, 
therefore, be interfered with each other. If the object is transparent, this allows sensitive 
detection of changes of optical path length through the object ( whether due to displacement of 
the object or changes in its Index of refraction. Similarly, two reconstructed wavefronts can 
be interfered with each other. This is most readily accomplished by super imposing on one 
hologram the wavefronts emanating from the object at two different times. The image recon- 
structed from this hologram contains a fringe pattern characteristic of any change in the op- 
tical path length through the object with respect to that at the time of the first recording. This 
latter technique was used in most of the experiments reported here* 

We have concentrated on the detection of optical path length changes which are due solely 
to refractive Index variation* The use of holographic interferometry to measure the refrac- 
tive index of homogeneous media has been explored, and an investigation of its application to 
heat and mass transfer and related fields has been initiated, 

o,2* BASIC APPARATUS AND PROCEDURE 

The apparatus of figure 14 can be utilized for interferometry in real time or for inter* 
ferometry by dual exposure* The procedure for real-time interferometry Is as follows: The 
photographic plate is exposed simultaneously to the reference beam and to the object beam, 
which has passed through the undisturbed transparent test object* When the developed photo- 
graphic plate is replaced in its original position,, the wavefront which emanated from the un- 
disturbed object can be reconstructed by illuminating the film with the reference beam. If 
the object beam is also turned on, it wall interfere with the reconstructed object beam- hence* 
if the test object is disturbed by displacing it or changing Its refractive index, fringe patterns 
are observed. In this manner f real-time differential Interferometry can be performed. 

Differential Interferometry by dual exposure is a similar process* In this case, two holo- 
grams are produced on the same photographic plate— one of the undisturbed test object and 
one of the disturbed test object. When this plate is developed and illuminated by the reference 
beam, the two wavefronts are simultaneously reconstructed and the resulting interference pat- 
tern is readily observed* 

In certain phases of the investigation, holographic interferometry was performed with a 
phase structure, such as a diffusing screen* placed between the collimating leas and the ob- 
ject. This gave rise to fringe localization phenomena of great Interest and possible applica- 
tion, which are described below* 

Experiments involving the measurement of refractive index and application to the study 
of natural convection heat transfer are described in the following sections* 
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5.3. MEASUREMENT OF INDEX OF REFRACTION OF HOMOGENEOUS MEDIA 

Conventional interferometric techniques lor determining the refractive index of transpar- 
ent objects entail considerable experimental difficulty. Generally, these techniques involve 
the comparison of optical paths In the two legs of a two-beam Interferometer. For the case of 
liquid or gaseous substances, each leg must include a transparent container — roe the substance 
of unknown refractive Index in one leg and for the standard substance to which it is to be com- 
pared in the other, U is essential that the two paths be Identical in the absence of the unknown 
and standard fluids. This requires precision alignment and high optical quality of all compo- 
nents, Including the containers. Once these conditions are satisfied, the resulting interference 
pattern provides an accurate measurement of the difference in optical path length between the 
two beams due to the presence of the test and standard llulds. The refractive index can then 
be computed. 

The use of holographic interferometry preserves the accuracy and sensitivity of the con- 
ventional techniques, but substantially alleviates much of their experimental complexity. This 
is true because the instrument is effectively a one-path interferometer. 

The holographic apparatus used in this experiment is the same as that of figure 14, ex- 
cept the collimating lens is removed from the object beam. Hence, the hologram Is the re- 
corded interference pattern between a spherical object beam and a plane reference beam. The 
test object in this case is a transparent container for the liquid whose refractive index Is to be 
determined. This container need not be or high optical quality; in fact, the one used in this ex- 
periment was made of plexlgLas. 

5,3,1, ANALYSIS, Let O denote the wavefront of the object beam which emanates from a 
monochromatic point source at fx^, y^, z^). Let r denote the plane wavefront of the reference 
beam which Impinges on the hologram at an angle 9 , The interference of these wavefronts re- 
sults in an Intensity distribution at the hologram plane given by 

f = IO + rj 2 - \Ol 2 + |rl a #0*r +Qr* 
where an asterisk denotes a complex conjugate. 

This Intensity is recorded on a photographic plate which is then developed such that its 
amplitude transmission at any point is proportional to I. When this hologram Is Illuminated 
with a monochromatic beam of wavefront C T a wavefront H. therefore, emerges from the holo- 
gram and is represented by 

H = C(fO| a + Jr ] 2 + 0*r + Or*j 

The first two terms within the parentheses represent the zero- order images and are of no in- 
terest in the present applications. The last two terms represent the first-order images. At- 
tention Is focused on the virtual Images Gr*C, 
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Let the coordinate system be located such that the hologram Is in the (*, y, O) plane. For 
simplicity, it is assumed that reference and signal -beam intensities are constant across the 
hologram, and also that the plane wave G used for reconstruction is identical to r. The re- 
sulting virtual image, which is a replica oi the object wavefront, is represented by 


where K is a constant amplitude term and where the phases of the reference and reconstructing 
beams cancel each other, 

Now suppose that a plane, parallel, and homogeneous transparent object such as a flat 
glass plate or a rectangular vessel containing a fluid is placed in the object beam. If the ob- 
ject is perpendicular to the beam, the point source will appear to have been displaced in the 
axial direction (see fig, 15), The resulting wavefront appears to emanate from a point source 
at (xy yj, z 2 ) and is represented by 


The interference of this wave with the reconstruction of the original object beam produces 
a Fresnel zone pattern of concentric circular fringes. When a first-order analysis is per- 


where A is the illumination wavelength* 

Since 2 ^ can be measured, the above relation can be used to determine Az - 2j - It is 
known that 


where n is the index of refraction of the object and t is its thickness; hence, if either of these 
quantities is known, the other can be determined readily. 

5.3.2. EXPERIMENTAL RESULTS, Experiments were conducted to verify the analysis 
presented in the previous section* Typical results for both dual -exposure and real-time in- 
terferometry are shown in figures IS through 18. The object was a transparent rectangular 
vessel whose parallel faces were normal to the object beam. In some cases, it was empty* 
Figure 18 is a photograph of the reconstructed wavefront from a dual -exposure hologram, the 
first exposure was made with no vessel in the object beam and the second exposure with the 
vessel in place. The expected Fresnel zone pattern is clearly evident. 




formed with and y 1 taken to be zero, the radius of the mth fringe in the observation plane 
is found to be 



Az - nt 
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FIGURE 16. DUAL -EXPOSURE INTER FERQGR AM FOR 
INDEX- OF REFRACTION MEASUREMENT 




(b) Recording of Hologram with Ho Object* 
During Reconstruction,. Object Inserted, 



(c) Recording ol Hologram with Object in Place* 
During Reconstruction* Object Removed* 


FIGURE 17. REAL-TIME INTERFERQGRAM3 FOR INDEX - 
OF -RE FR AC TION M E ASUR EM EN T 
45 
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fa) Effect of Emulsion Distortion 





(b) Recording of Hologram with Object in Place* 
During Reconstruction, Object Removed* 



(c) Object Illuminating Source Moved to Null Position 

FIGURE IB. REAL-TIME INTERFEROMETRY FOR INDEX- 
OF-RE FRACTION MEASUREMENT— NULLING TECHNIQUE 
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Analogous results of real-time holography are shown in figure 17* Figure 17a shows the 
inter lerence of waves from the actual point source with those from Its own reconstruction. 
Although inaccurate positioning of the hologram and distortions of Us emulsion caused by the 
development process made it impossible to achieve a true aero fringe situation, a very low 
order of fringes was obtained* The vessel was then inserted into the object beam and the re- 
sulting fringe pattern is shown in figure 17b. The Inverse procedure was also carried out: a 
hologram was made with the vessel In the object beam. It was then removed during recon- 
struction. The resulting fringes are shown in figure 17c. As expected, the results of all three 
holographic procedures are equivalent* 

A different technique which allows direct measurement ol optical path length changes was 
also tried, A hologram was constructed with the water-filled vessel in the object beam. Alter 
appropriate development, the hologram was returned to Its original position. The resulting 
course fringes are shown in figure 10a. The water was then removed from the vessel and the 
resulting fringes at the observation plane are shown in figure 10b. The point source was then 
axially displaced by means of a micrometer-driven platform until the fringes were nulled. 
Although we could not obtain the one fringe situation in figure 18a, a rather low-order fringe 
pattern was reached (fig, IBcJ. The displacement gave a direct measurement of &z t from 
which the refractive index of the water could readily be calculated. With our relatively crude 
setup, we were able to achieve 5% accuracy with relative ease. A more sophisticated appa- 
ratus could no doubt greatly increase this accuracy* 

5,4, EXPERIMENTS WITH NONHOMOGENEOUS MEDIA 

5*4.1* DIFFERENTIAL INTERFEROMETRY WITH PLANE WAVES, To illustrate the 
use of holographic interferometry of nonhomogeneous media, we examined the density' field of 
air surrounding heated objects. Two soldering irons, about 5 -in, apart* were placed In the ob- 
ject beam of the apparatus shown In figure 14, An exposure was made while the irons were in 
thermal equilibrium with the ambient* The iron closest to the hologram was then allowed to 
heat for 5 min and a second exposure was made on the same film plate* The fringe pattern 
observed upon reconstruction of this dual-exposure hologram is shown in figure 10a. These 
fringes are non localized: i*e*j in focus everywhere. The photograph was made by simply ex- 
posing a sheet of Polaroid film to the beam emanating from the hologram — no camera was 
necessary* Similar results were obtained by the real-time procedure. The fringe patterns 
observed In this case are entirely analogous to those produced by classical interferometers 
such as the Mach-Zehnder* 

A second experiment was carried out in the same manner, except this time both irons 
were heated simultaneously. The resulting nonlocallzed fringes are shown in figure 19b. 
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(a) One Iron Heated 



fb) Both Irons Heated 

FIGURE 19, INTER FERQGFE AMS OF NATURAL CONVECTION 
IN AIR — PLANE WAVE INTERFEROMETRY 
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5.4,2* DIFFERENTIAL INTERFEROMETRY WITH DIFFUSED LIGHT. A second group 
of experiments, like those In section 5*4,1, was carried out, but this time a diffusing screen 
was placed in the object beam such that the light passing through the object was diffuse. The 
Tringe patterns observed upon reconstruction of the various dual-exposure holograms are 
Shown in figures 20,. 2 1 r and 22. These fringes, unlike those discussed in the previous section, 
are localized at the position of the heated irons. When viewing these holograms, we found it 
readily apparent that fringe patterns shift as the viewing angle is changed. This localization 
phenomenon is quite obvious in the sequence of photographs shown. 

Figure 2 Da shows the fringes obtained when only the Iron farthest from the hologram was 
heated. Figure 20b shows the analogous fringes when only the iron closest to the hologram 
was heated* Figures £la and 21b are two different photographs of the same hologram. In 
this case, both irons were simultaneously heated. It is obvious that the density field around 
each iron can be individually visualized— a phenomenon which has no analog in classical in- 
terferometry. Finally, figures 22a and 22b are photographs of the reconstruction of a triple- 
exposure hologram* The first exposure was made with both irons in equilibrium with the am- 
bient. A second exposure was made with only the iron closest to the hologram heated, and 
later a third exposure was made with only the Iron farthest from the hologram heated. Again 
It is obvious that we can focus on the individual fringe patterns, 

5,5, EVALUATION 

The experimental results Indicate that holographic interferometry is a potentially impor- 
tant tool for the measurement of the retractive index of phenomena which manifest themselves 
as inhomogene It ies of transparent materials* Examples of such phenomena are aerodynamic 
flows, convective heat transfer* mass diffusion, and stress analysis. One of the early dramatic 
examples of this kind of Interferometry applied to transparent media is, of course, the shock- 
wave patterns of a bullet In flight demonstrated by Brooks, He! linger, and Wuerker (in Applied 
Physics Letters , Vol* 7 , 1955, p* 92). 

Although most of these measurements can be made by classical interferometric methods, 
the relative ease of setting up and adjusting a holographic interferometer may make it a pref- 
erable instrument for many applications, A particularly attractive feature of this method is 
the elimination of the need for compensating chambers and for test sections of high optical 
quality in heat transfer and aerodynamic experiments. 

The Interesting prospect of interferometry with three-dimensional density fields is raised 
by the fringe localization effects observed when diffuse object beams are used. Evaluation of 
this prospect* however , awaits further investigation and understanding of the phenomenon. 


(a) Iron Nearest Hologram Heated 



(tij Iron Farthest from Hologram Heated 


FIGURE 20, INTER FEROGR AMS OF NATURAL CONVECTION 
IN AIR -DIF FUSE ILLUMINATION INTERFEROMETRY 



(a) Interfere gram of Nearest Iron to Hologram 



(b) Merlerogram ol Farthest Iron from Hologram 

FIGURE 2L INTERFEROGRAMS OF NATURAL CONVECTION 
IN AIR -DIFFUSE ILLUMINATION INTERFEROMETRY— IRONS 
HEATED SIMULTANEOUSLY 
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(a) Interf program of Nearest Iron to Hologram 



(b) Interior ogr am of Farthest Iron Irom Hologram 

FIGURE 22. INTERFERQGRAMS OF NATURAL CONVECTION 
IN Am- m FT USE ILLUMINATION INTERFEROMETRY— IRONS 
HEATED SEPARATELY 
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6 

PROPAGATION THROUGH A TURBULENT MEDIUM 

6.1. INTRODUCTION 

It has been repeatedly observed in certain propagation experiments (see fig. 23) at optical 
frequencies that the first-order distribution of the received wavefront is log-normal. We 
would like to show the plausibility of this fact and to investigate a method of approximate solu- 
tion of the wave equation that is suggested by this experimental result. 

The log-normal distribution is, of course, directly predicted by the geometric optics and 
Rytov approximate solutions of the wave equation* (The latter method has recently fallen into 
some disrepute [16, 17) and, therefore, cannot be considered a reliable predictor.) On the 
other hand, the Born approximation, of course, amounts to taking only the first two terms of 
the Neumann series solution. Ignoring convergence questions, the resulting distributions of 
the scattered wave (as given by the Neumann series) will be quite complicated as the sum of 
all order products, roughly speaking, of the refractive index field are involved. 

We shall review the basic way the log-normal law occurs, and this will suggest a method 
of approximate solution of the wave equation. The method will involve certain iterations of 
the Born approximation. In order to end up with a "closed” form for the approximate solution, 
we will have to make certain assumptions about the nature of the scattering and the commuta- 
tively of certain operations. 

We do not believe the method of solution is new, although we are not aware of similar ap- 
proaches for optical waves in the turbulent atmosphere. The closest work appears to be that 
of Feinstein [18] who calculated the covariance function of such a field. The method appears' 
to be in the vein suggested by Bellman and Koloba [19] and, although for "discrete" scatterers, 
Feynman [20]; in particular, the latter employs formally an "operator calculus" as we do here. 

6.2. THE LOG-NORMAL DISTRIBUTION 

The motivation for the present discussion is the need for an adequate characterization of 
an optical wave propagating a relatively long distance through a turbulent medium represent- 
ing the atmosphere in order to eventually find the statistics of the scattered field. In turn, 
the statistics of the scattered field are required in order to find optimum procedures for 
transmission, reception, and processing of information transmitted via the atmosphere by 
waves in the visible spectrum. Of course, it is generally true that more second-order statis- 
tics are required to determine optimum procedures. 

We emphasize that our discussion is heuristic and, again, we are led to it by the empir- 
ical observation of log-normal statistics rather infrequently; the log-normal probability law 
occurs for a very natural reason. Suppose a random phenomenon is described by a family of 
random variables {E n } such that, for every n, 
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. FIGURE 23. TYPICAL GEOMETRY FOR OPTICAL COMMUNICATION CHANNEL 



(32) 


E ,, = E + 
n+1 n 


W E n 


where the {^} is a family of random variables. The meaning of equation 32 is clear: the 

n + 1th random variable is the sum of the nth random variable plus a small perturbation which 

depends on the nth random variable; specifically, the dependence is direct proportionality. The 

fact that many natural phenomena can be modeled this way is not surprising, and the application 

to the progression of a wavefront from a plane P at z = z through a turbulent medium to a 

n n 

plane at z = is already roughly indicated. 

We assume that ^ is small in some appropriate sense, e.g., | « 1 with large proba- 

bility. Equation 32 can be written 


and hence. 



(33) 


log 


“n+1 


n+1 

£=0 


(34) 


n+1 

^J2+l is (tends, as n — co, to) a Gaussian, or normal random variable, then (by definition) 

log (E^^/Eq) is (tends, as n — co, to) a log-normal random variable. 

There are a variety of conditions that allow a normal approximation to this summation: 

(1) If the are normal, then so is every partial sum. 

(2) If the {£ n } are independent and identically distributed with mean m and variance <x 1 2 , 

n 

then as n — co, the partial sum converges (in distribution) to a normal random 

1 2 

variable with mean nm and variance no (see the Lindeberg-Levy theorem in ref. 21, 
p. 215). 

(3) If the {i;^} are independent but not identically distributed, there is a central limit 
theorem in agreement with Lindeberg [22, p, 236] . 

(4) There are central limit theorems for partial sums of dependent ^variables [23, ch. 8], 

The above discussion can be easily extended to continuously varying perturbations 
x now plays the role of a continuous index and £(x) is a random process. Now, for a small 
change A in x, we have 


E 


x+A 


E + A£ 
x 


and, as n — oo and A — 0, we find 
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(35) 


log 




i dx 
x 


It is interesting to note that equation 32 implies a specific way of accounting for multiple 

scattering. Consider Eg = (1 + ^ E 1 = ^ + + ^ E 0 = (^ + + ^ + ^1^P E 0* This scat " 

tering is shown graphically in figure 24: we see that only "forward” scattering is accounted 
for. This is roughly the kind of approximation that we want for the scattering of optical waves 
by atmospheric turbulence since it is known that, approximately, only forward scattering in 
small angles is possible. 

It should be pointed out that there are situations, different from the geometry of figure 23, 
when a log-normal law should not be expected. For example, consider the situation shown in 
figure 25: If the receiving aperture is far from the boundary of the turbulent region, relative 
to the correlation distance of wavefront emerging from the turbulent region, then the field at 
the receiving aperture is approximately a summation of uncorrelated contributions. Though 
this is not sufficient, we might hazard that the field at the aperture is approximately normal — 
not log-normal. 

Generally there would be the intermediate situation where contributions to the field at a 
given point represent a sum of the contribution propagated along a ray from source to the 
point — which can be log-normal — and a sum of the contributions scattered from other rays — 
which could be normal. 


6.3. APPROXIMATE SOLUTION OF WAVE EQUATION 

As is well know [24] , the inhomogeneous scalar wave equation in a medium with varying 
2 

index of refraction n = 1 - q is 

V 2 E + k 2 E = E q +k 2 qE (36) 

provided the index changes very little over a wavelength. In (36) we have assumed, as usual, 
a sinusoidal time variation for the field; Eq is the scalar field of the source. If satisfies 
homogeneous boundary conditions, the radiation condition (for unbounded media), and the equa- 
tion 

V 2 E 6 +k 2 E 6 =6 (37) 

then, (36) may be written as an integral equation: 

E =E q * (e q + k 2 qE^ s AE Q + k 2 A(qE) (38) 

We have introduced the operator notation A: Af h E^ * f (A is a linear operator). The identity 

operator will be denoted I: K = f . The Born approximation Eg to E is the first two terms of 

the Neumann series [25] which is generated by simple iteration: 
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FIGURE 24. THE SCATTERING MECHANISM, IMPLIED BY EQUATION 32, THAT IMPLIES A 

_ n 

A LOG-NORMAL DISTRIBUTION FOR E^/Eq AS n - oo IF^~~| ^ CONVERGES (IN DISTRIBUTION) 

1 

TO A NORMAL RANDOM VARIABLE AS n - oo 
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FIGURE 25. GEOMETRY WHERE NORMAL, VICE LOG-NORMAL, DISTRIBUTIONS MAY OCCUR 
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(39) 


E b = AE fl +k 2 A(qAE 0 ) 


Referring to figure 26, we wish to calculate the field after it has propagated a distance L 
into the turbulence. We allow L to be a distance too large for the Born approximation (Eq. 39) 
to be good. We subdivide the region (0 ^ z ^ L) into layers A thick as shown (z^^ - z^ = A, 
all n) such that — the Born approximation is good, and the free -space diffraction is approxi- 
mately 


AE 

n 


i(z-z ) 
n ' 

e E 


n 


(40) 


With equation 40 and the fact that scattering will only be into small angles, we can make the 
simplification 


ik(z 1 -z ) fff iktf> ik(z -z ) 

A(qAE n ) « e Jjj fafe - z ) q(r)E n (x > y)dxdydz “ 6 A l ( VlV (41) 

V, n 


where 


, m \ 


n+1 - X)2 + ^n+l 


2{ Vi - z) 


,21 

- y) 


the so-called sagittal approximates. Then equation 39 becomes 


k „*i * ‘““t 1 *‘\vK * A ‘ <Vl \ 


(42) 


Equation 42 is the generalization of (32) for the propagation of an optical wave through a 

turbulent medium of described type; exp k B is, of course, a linear operator defined by the ex- 
2 

pansion 1 + k B + „ . . . We index q because it depends upon the particular sample function 
segment of q present in the slab {r: z n ^ z ^ We iterate the relation (32) and find 


E « e 
n+l 


iknA 


n 

TJ e 

£=0 


k 2 A, 


2 ^+ 1 ^ 


E, 


(43) 


The product of operators appearing in (43) is generally ordered. 

However, in the case of homogeneous turbulence, that is, turbulence whose statistical 
properties are invariant to spatial location, any particular segment of q appearing in the 
slab {r: z. izi z^} is i ust as likely to appear in the slab {r: z. ^ z £ j £ i. There- 

fore, the statistical properties of E ^ should be invariant to permutations of order in equa- 
tion 43. We therefore assume that, insofar as we are interested in statistical properties of 
E p the operators may be assumed commutative. But if this is true, then it is easy to show 
that 
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FIGURE 26. THE GEOMETRY AND SUBDIVISION FOR THE CONVOLUTION OF THE 

SCATTERED WAVE 



where 


» kVv r > 

// e = exp 

£=0 


E a i ( Vi 


k' 

L £=0 


= exp |k^A 


[“Vv-i 


jmj 


Thus, finally, 


„ ^iknA 
E n+1 * e e 


ik‘ 


' A i«w> 


E, 


(44) 


6.4. APPLICATION 

We now return to the geometry of figure 26 and apply the above approximation method. 
We. divide the z-segment {0 ^ z ^ L} into N = L/ A "slabs” where A is chosen so that (1) the 
Born approximation is good, (2) free-space diffraction is small, and (3) the segments of q 
are approximately independent. Requirement 2 will be satisfied if, roughly, 


where & is the smallest significant dimension of the turbulence. Requirement 3 will be satis- 
fied if A » L q where L Q is the largest significant dimension of the turbulence. We must, 
therefore, have 



» 1 


This is presumably true in the atmosphere for optical wavelengths. 
Then, we have 


where 




(45) 


»*■ 

£=1 

and A a is defined by equation 41. Equation 45 is the result that we have sought. 
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6.5. FURTHER COMMENTS 

The result (45) states, among other things, that the turbulence is equivalent to a strange 

2 

turbulence distributed over a correspondingly shorter distance. It is to be noted that k A A (q*E Q ) 
•need not be small (in some sense) with unity in order that (45) hold. For example, in the ray 
optics limit, 

ijJ q*(z)dz 

A A (q*) - e A (46) 

the phase need not be small with respect to a radian in order that this approximation apply. 

If the region {z: 0 ^ z ^ L} is composed of subregions of homogeneous turbulence, of a 
different statistical nature in each region, then, of course, we can apply our approximate method 
of each subregion. Obviously, some of the subregions can be free of turbulence. (We must 
keep in mind that we have ignored free- space diffraction: this limits the distance in a 
well known way on beam dimension and wavelength. For example, for a uniform circular beam 
area A, we must have « A/x.) 

It should also be noted that result (45) does not necessarily predict log-normal statistics. 
Generally, the statistics are quite complicated. In the ray optics limit, we generally expect 
log-normal for the reasons originally outlined. 


7 

THE LOG-NORMAL RANDOM PROCESS 


7.1. DEFINITIONS 

Both the experimental results on the first-order distribution only (apparently) and the 
Rytov method within its domain of validity show that the so-called log-normal distribution oc- 
curs in optical propagation through the turbulent atmosphere, although the exact conditions of 
its appearance are unknown. This process is, therefore, worthy of study, and is analyzed here 
in relation to the definitions listed. 

(1) Let (X fc , t e T), X fc - + iS t be a complex Gaussian process; then the process where 

{Y t , t £ T} and 

x f 

Y t s e , t € T (47) 


is called a log-normal process. 

Definition (1) is a special case of the more general nonlinear transformation g(x fc ) of {X t , 
t e T}, where g is a real-valued measurable function of a real variable. 

Equation 47 is such that the ensemble averages of combinations of the values of Y^. can be 

3 

evaluated from the well-known forms of the characteristic functions of a Gaussian process. 
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For the optical propagation problem, we wish T to be of greater than one dimension; then 
'’process” is replaced by "field” in definition (1). 

(2) We recall that if (x n ) is an aggregate of real Gaussian random variables, then its joint 
characteristic function is 

<f> x (Uj, . . . , u n ) s E[e^ u * x ^] = exp ^iu^m - 

where the row vector u* = (u^, • . • , u N )isthe transpose of the column vector u, m, is the column 
vector of means, in* * (m^ . . . , n^O), m n = E(x n ), and A is an N x N covariance matrix of 
element A mn = E{[x m - E ( x m )H x n “ E ( x n )]}* Their joint density function is 


N 


P (x y ^ = 1 1 

v > ♦ * * ? v v*in> a * * A */nv _ i — — r ex P 


T" 1 I A | (x n - m )(x n - m ) 

/ , mn v mO m v nO n 


x m io ,- ‘ M0/ 


ttat 


(48) 


where I A | is the determinant of A and | A (is the determinant of the cofactor of A * The 

mn mn 

cofactor matrix A_^ of A is obtained from A by omitting the mth row and nth column of A. 
mn mn J b 

The statements of definition (2) carry over immediately to a random process (x^_, t e T) 

with mean- value function m(t) s E(x.) and centered covariance function r(t, s) = E[(x f - m,)(x 

t t t s 

m s )]. 

(3) For any complex process (z^, t e T), we define the following expectations: The mean- 
value function m (t) = E(zJ; the covariance functions R (t t , t 0 ) s E z, z, and R (t 1 , t 9 ) = 

( \ Z L Z i 6 1 / ZZ J. cx 

S\} 

We have R z (tj, t 2 ) = and t^) = R- zz (t 2 , tj). 

If z„ s Re (z.) and z T = Im(Z), then we define the covariance function 
R t 1 \ 

R z z (t l’ V S E ( Z R Z I 

R i \ tj t gy 


and have immediately 


R z R z ][ (t l , . t 2 ) " R z^ r ^2» V 


Then, 


R A- V ■ V + V»’ V * S) - l Vl (t ’ S) 


and 
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R zz^l’ V " R z r ^1’ V ' V + ^z^’ s) * lRz R z i (t ’ ^ 

We shall also use the centered covariance functions 

- = R z (t r t 2 ) - m z (ti ) m z (t 2 ) 


* z h> V s E 


\ - m z ( V 


\ - ra z ( V 


and 


r zz (t l’ fc 2 } B 


\ - m z (t l } 


\ " m z ( V 


= R zz (t l’ t 2 ) - m z {t l )m z (t 2 ) 
Lemma: If [Y t s exp + iS t ), t e T] is a log-normal process and (t n e T, 1 g 


(a , 1 ^ n ^ N) are complex scalars, then 


where 


' a l a N\ iu t m-iu t Au 

E|Y. . . . Y, =e 
1 N/ 


u - (u l; . . . , ^2N) 

-icy - , 1 s j ^ N 
3 

U. H*' 

- 1 |o^ N + 1 ^ j ^ 2N 

ni = (m^, * * ♦ 1^2 


|mp(t.), 1 - j -N 


m. =- 
3 


m q (t.), N + 1 * j ^ 2N 


and A is the 2N x 2N matrix of elements 


X = E[(? - m )(£ - m )] 

mn LK m m v n n /J 


€ m s ' 
m 


m m ), 1^0 

|S(t ), N + 1 ^ m ^ 2N 


Proof: We have 


/“ 1 “n\ 

+ • * * +i(_ia N )(3 t i “i S t 1 + ’ * • + i Vt N 

h-* 

• 

II 

H 

e 1 N -e 


n ^ N) and 


(49) 
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which is the joint, characteristic function of the jointly normal variables f/3. , . . . , j3, , S f , 
\ \ h *N 

* • • > s t ) evaluated at the point (-ic^, . . . , -ia N , . , a N ). Equation 49-then follows 

directly from definition (2). 

Theorem 1: If [Y^ = exp {/^ + iS^), t e T] is a log-normal process and {t^ e T, 1 ^ n ^ N} 
and a = (a j, . . . , a^) is a row vector of N complex scalers, then 


E 


v\ ‘ "\) = exp [“H + + ~ *** + + -V)j 




where is the row vector of elements (t n ) = E[|B(t n )], nig is the row vector of elements 


m g(t n ) = E[S(t n )j, A ^ is the N x N matrix of elements 


r ... (t , t ) s Ei 
4>v m’ n 1 


h - E 

m \ 


m 




t 


where and are either por S. 

Proof: We note that the vectors and matrix of the lemma can be partitioned 


u* =5 (-ia z | £a L ), m L = (m^ |m g ) 




and 


A__ 

A „ 

JSJS 

j3S 

A S/3 

^ss 


with definitions given below in the statement of the corollary . Since matrixes with partitioned 
matrix elements obey the same multiplication rule as matrixes with scalar elements (order 
must be preserved, of course [26]), the assertion now follows. 

Corollary 1: If (Y fc = jS^ + iS^, t c T) is a log-normal process, a is complex scalar, and 
t e T, then 


E (^“) =exp + i“m g (t) +i[o- 2 CT^(t) -i-2o! 2 cr^(t)cr s (t)p^ g (t, t) - o: 2 cr s (fc)]] 

9 9 

where Ppg(t> s) s r^g(t, s)/a^(t)a g (s), o* s r^(t, t), and n g e r g (t, t); in particular, the mean- 
value function 


m y (t) S E(Y t ) = exp jm^ + im^ + f o 2 ^ + 2*^ p^(t, t) - j 

Proof: In theorem 1, the vectors a , m^, and m g reduce to scalars a ^m^(t), and m g (t), re- 
spectively, and the matrixes A^, A gg , and A^ g reduce to scalars (t, t), r gg(t, t), and 


r^g(t, t), respectively. 
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Corollary 2: If [Y^ = exp 03^ + iS^), t e T] is a log-normal process; a 2 are real num- 


“l“2 


bers; and t^, <~ T, then E jY^. Y t J is given by theorem 1 with 

“ t S = “ 1 m / 3 (t 1 )+ “ 2 m | 3 (t 2 ) 


am s = a 1 m s (t 1 ) + “ 2 m S^V 

“ A /3^“ = a l CT |3^1^ + 2 “l a 2 ff (3^1^ CT /3^2^ P i3^1’ V + a 2°VV 
a A SS a = a l a sh ) + 2 a l a 2 a S^ ff S^ p S^l’ V + a 2°W 
01 A j 3 S a = “l r j 3 S (t r V + i a 2 V BS^l’ V +a 2 r 8S <<i "2’ V 


a A Sp° L '~ a l X Sp( t l > V + 2a l a 2 T S^ t l’ V + a 2 r Sj3^2’ V 

In particular, the covariance function R^-y is obtained from these relations with = a 2 s 1. 

Corollary 3: If [Y^. s exp 03 ^ + iS^.), t e T] is a log-normal process, then its covariance 
function is 


Mr t 2 ) s E ( Y t 1 Y t 2 J = exp pj + m (3 2 + i \ m S 1 - m s. 


+ i R +2R +R - R -2R +R 
2 \ ^11 p 12 p 22 / \ S 11 S 12 fa 22 . 


+2i \V n " v 2a 


Proof: The proof is straightforward but slightly involved* First of all, 


EY tM =* 


1 2j ^,8^2 


(-i,-i,l,-l) = exp 


where 


U* B ^u* |u 2 ) = [(-i - i) I (1 - 1)] 

= ( m fp 2 ) (“sp, 
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and 


where, e.g., 


a./A 


i A, 


SJ3 


(SS 1 


A, 


R, 

A -I % ^ 

0 l Bp B 

P 21 P 22 , 


With these partitionings, we can write 


^Au = S*A/i + u^ s u 2 + ^A^ + u*AgU 2 


and 


— *■! ***■ ~ ~ 

U m=u l m j3 +u 2 m S 

By direct calculation, we find, e.g., 

ui.Au = - R +2R +R 

1 P 1 V *11 ^12 h2j 

where we have used the fact that R = R for a real process (|3. ). Further, 

p 12 P 21 1 

"^S”2 -IVl ' -‘fsjj, - E S 0 j 

where we have used the fact thar R flQ = R„ fl for real processes. With similar calculations, 

^ S 12 S ^21 

the assertion results. 


7.2. STATIONARY PROCESSES 

We assume, in this section, that (p^., t e T) and (Sj., t e T) are wide-sense stationary, and 
we list the following definitions: 

(4) The real random process (x^, t e T) is wide-sense stationary, if e(x^ < co and R x (t + 
s, s) s E(x t+g x g ) does not depend on S. 

If the real wide-sense stationary random process {x t , t e T} is Gaussian and E(x t ) does 
not depend on t, then (x £ , t e T) is strictly stationary. 

(5) The real, wide-sense stationary processes {(3^, t e T} and {S t , t e T} are jointly 
wide-sense stationary, if R^t + s, s) = E{|3 t+g S g } does not depend on S. 
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We confuse the notation by denoting, e.g., R_(t + s, s) by R^(t) when t e T} is wide- 
sense stationary. If {|3 t , t e T} and {S t> t e T} are jointly wide-sense stationary, then V« 

- V' t% 

(6) A complex random process (X^., t e T) is wide-sense stationary, if E |jx t | 2 j < co and 
R (t + s, s) e E{X. X } does not depend on s. 

From equation 48 above, we see that (X fc s j3 fc + iS t ) is wide-sense stationary, if {fy, t e T} 
and (Sj_, t e T) are jointly wide-sense stationary. 

If a complex process is wide-sense stationary and Gaussian with, e.g., zero mean value 
function and R^t + s, s) s E{X t+g X g } s 0, then it is strictly stationary. 

We shall assume where the process (fy., t e T) and (S fc) t e’T) are jointly wide-sense sta- 
tionary, that their mean value functions (t) and m g (t) do not depend on t. 

Theorem 2: If {Y fc s exp (fy + iS fc ), t e T} where {fy, t € T} and {S t ,t e T} are jointly wide- 
sense stationary Gaussian (JWSSG) processes with constant mean value functions, then 

m y = exp jm^ + im g + 1/2 [a 2 + 2 g^j (0) - or J} (50) 

R y (t) = exp |2m^ + ff 2 [l + (1)3 - <r- g [l - p g (t)]j (51) 

and 

R yy (t) = R y (t) exp 2i[m g + r^(t) + r^(t) + 2r^(0)] (52) 

Proof: Equation 50 follows trivially from the equation of corollary 1. Equation 51 follows 
easily from the equations of corollary 2. Equation 52 follows trivially from corollary 3. 

Theorem 2 implies that if {/3 fc , t e T} and {S t , t e T} are JWSSG processes with constant 
mean value functions, then [Y fc s exp (fy + iS fc ), t e T] is a WSS log-normal (complex) process 
with constant mean value function. 

Theorem 3: If [Y t s exp (£ fc + iS^.), t e T] where (f3 fc , t e T) and {S t , t € T} are JWSSG 
processes with zero means, then the covariance functions R y and R^ y completely determine 
the (Y £ ,’t e T) log-normal process. 

Proof: It is sufficient to show that R and R determine R„ and R , for then the Gaussian 

y yy p ° 

processes {/3 t , t e T} and {S^., t e T} will be completely determined (recall the lemma) and 
hence then also (y^, t e T). 

From theorem 2 we have 

log R y (t) = ff|[ 1 + p^(t)] - ff g [l - P s (t)] (53) 

and 
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( 54 ) 


lOg = CT|[l + P^t)] • 0^1 + P s (t)) 

Subtracting equation 54 from equation 53, 

2 , E v« 

^p s (t) = 1/2 log 

yy w 

thus, (s^., t e T) is completely specified* 

We can then, e.g., add equations 53 and 54 to get 

1/2 lOg Ry(t)Ryy(t) = 0-|[l + P^(t)] ~0 2 s 


(55) 


or, 


ap + p/t)] = log R y (t) 

Hence, it is easy to see that 

2 R .,« 

Thus, t e T} is completely determined. 

Lemma: If t e T) and (S^., t e T) are WSSG processes with constant mean value func- 
tions, then 


E e 


WO. 


= exp|2m^ + CF^[1 + p/t)]} 


(56) 


and 


i{ W iS sl 


ex P {-o-gf 1 - P s (t)]| 


(57) 


Proof: To prove equation 56 we, of course/ note that E[exp (j3 + j3 )] is the joint charac- 

ters s 

teristic function of the Gaussian random variable (3. and /3 , evaluated at (~i, ~i). From the 

t+s s 

lemma, we have the asserted result. Similarly, E[exp (iS, - iS )] is the joint characteristic 

t~rS S 

function of the Gaussian random variables and S g , evaluated at (1, -1). From the lemma, 
we have the asserted result (eq. 57). 

Theorem 4: The log-normal random variables exp 03^ + /3 g ) and exp i(Sj. +s - S g ) formed 
from the JWSSG constant mean processes t e T) and (S^, teT) are uncorrelated. 

Proof: We note that we can write 

R y (t) = E[exp G3 t+£ . + P g )-exp i(S t+g - S g )J 
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Comparing the form of R y asserted by theorem 2 with the asserted forms of the lemma, we 
have 

E[exp (0 t+s + f3 g )-exp i(S fc+s - S g )] = E[exp ® t+s + /3 S )]-E[exp i(S fc+s - S g )] 

Theorem 4 has been asserted with the incorrect proof that asserts the two random varia- 
bles (p. + |3 ) and i(S. - S ) are uncorrelated. For example, Cernov [27] specifically 

^l+S s t+s s 

asserts E(0 t2 S tl ) = E(f3 tl S t2 ) which is not generally true. The most that can be said without 
further assumption is that R^g(t) = Rg^(-t) (they are real, WSS processes). 


7.3. LIKELIHOOD FUNCTIONAL 

Of paramount importance in communication-type problems of detection (hypothesis testing 
and estimation, e.g., regression) is the likelihood ratio, or likelihood functional, in the context 
of random -processes. We now seek this functional for the log-normal random process. 

Real processes: zt s exp x fc . We begin by considering real log-normal processes {z fc = 
exp x fc } and the typical test between the two hypotheses: 

H l : z t = g(x t } 

H 0 : z t= g(y t ) 


where (x^., t e T) and (y fc , t e T) are Gaussian processes. We suspect that if g is, e.g., a mono- 
tone function of a real variable, the likelihood functional should exist provided the likelihood 
functional for the problem 

H r s' lz t =x t 

H o : g "S =y t 

exists: that is, if the Gaussian processes {x^, t € T} and (g^., teT) are equivalent. 

The basis for this expectation is the following. We take as the underlying probability 
spaces the function spaces: then g is a very simple transformation on these spaces, and the 
new measures seemingly can be obtained from the old just as with ordinary random variables. 


f z [g~ 1 (z)]^g" 1 (z) yg' 1 ^] 

yg-^Ag-lfc) 'yg^z)] 


=f 1 [g _1 ( Z )] 


We may motivate this further by the following argument. Under conditions of equivalence, the 
likelihood functional is calculable from the limit 




> Z tN^ 

» Z tN^ 


70 



To find this ratio for when z = g(x) and for H Q when z = g(y), we use well-known' elementary 
probability theory. The Jacobian reduces to a diagonal because of the simple nature of the X 
formation from x(y) to z: 



In any case, the Jacobian cancels out, giving 


W z) = 


f 

X 

£ ( z ti^ • * • 9 S 

f 

y 

S • ■ • > ( z tN^ 


It would appear that under certain conditions of g, the linear operation should be successful. 
Formally, 


,(z.) = lim f„„(z) =f 1 g 1 (zJ 
1 L N-*-co tJ 


^lt ^9 1 

Theorem 5: If = e + A^ = e are log-normal processes such that ((3^) and (f^) 


It 


are equivalent, with likelihood functional f^(z) ■ dP^ (z^dP^ (z), then the likelihood functional 


f A< z) 


dP Al (z) 

dP^J 


- IpOog z t ) 


Proof: We have the measurable space (ft, j3) associated with the process (x t ,t e T); on the 
measurable sets we have either a probability measure P Q or a probability measure P^, giving 
rise to two measure spaces (ft, ]3, Pq) and (ft, j3, P^), We assume that Pq and P^ are equivalent, 
i.e., mutually absolutely continuous. Then, by the Radon-Nikodyn theorem there exists a finite- 
valued Pq measurable function f such that, for any set B e /3 


P^B) = 


1 


f(w)dP 0 (w) 


(58) 


f is the Radon-Nikodyn derivative, frequently denoted dP^/dP^, or likelihood functional [28, 
theo. 31.B] . 

Now let g be a monotone, increasing, continuous, real function of a real variable and con- 
sider the new random process, defined by the relation y(t) h g[x(t)], with measurable space 
denoted (2, S). This relation induces a transformation T from O to 2: its nature is most 
easily seen if we take and 2 to be the function, or coordinate, spaces 5 and 2 [29]. For 
then T is g: since g is continuous and monotone, T must be measurable and one-to-one. Thus 
T induces two new measures, Qq = P^T * and Qj = P^T - *, giving rise to the two measure 
spaces (2, S, Q q ) and (2, S, Q^. 
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The measures Q q and are also equivalent, since P Q and Pj are [28, p. 165]. There- 
fore, there exists a Qp-measurable function h such that, for any s e S, 

Q 1 (S) = J h(o)dQ 0 (o) 


That is, 


•i 1 ' 1 ® - 1 


h(ff)dP 0 T _1 (?) 


But we can rewrite the right side as [28] 


| h[T<«)] 

m 1 /rt\ 


dP 0 (w) 


T A (S) 


thus, setting B = T _1 (S), we have that for an arbitrary s e S and hence B e 0 that 


?1<B) - f 


h[T(5)]dP 0 (£) 


But comparing (58) and (59), we have f(w) = h[T(w)]; that. is, the Radon-Nikodyn derivative 
(likelihood functional) for the (y^., t e T) process is 


Since this can be written 


h(a>) (T _1 co) 


h(x t ) =f‘^ 1 [g(x t )] 
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Appendix 

PHOTOCHROMIC MATERIALS IN HOLOGRAPHY 

A.l. INTRODUCTION 

During the past year an extensive program was designed to investigate the capabilities of 
photochromic materials as three-dimensional recording media for holographic data and to in- 
vestigate their utilization in real-time processing. These materials included F-centered ' 
potassium bromide (KBr) crystals, silver halide glass, and a variety of photochromic plastics. 
A significant amount of data has been obtained that is directly applicable to construction of 
holograms, including 

(1) Exposure sensitivity 

(2) Angular orientation sensitivity 

(3) Diffraction efficiency 

(4) Wavelength sensitivity 

(5) Reciprocity failure 

(6) Temperature effects 

This section presents a large part of these results and also contains information on our expei 
iments with multiple -image storage and multicolor wavefront reconstruction. 

A.2. EXPOSURE SENSITIVITY 

One of the properties of photochromic materials that makes them potentially useful is the 
occurrence of a change in optical transmission during exposure to a source of radiation of the 
proper wavelength. Because this change in transmission is concurrent with the exposure 
process, no developing is required as is usually the case in photographic emulsions. In addi- 
tion, this change can be reversed by reilluminating the photochromic material with a different 
wavelength. 

It is important to determine the amount of energy that is required to cause a certain 
change in the optical transmission. These changes, or exposure characteristics, can be pre- 
sented in several ways: amplitude transmittance vs. exposure, density vs. exposure, and den- 
sity vs. logarithmic exposure. This report presents plots of amplitude transmittance vs. ex- 
posure which can be used to determine the sensitivity (time of exposure), the region of linear 
response, and the operating point which obtains maximum diffraction. 

The exposure sensitivity depends on many factors, including (1) type of photochromic ma- 
terial, (2) wavelength of illumination, (3) temperature, and (4) intensity of illumination (reci- 
procity). Amplitude transmittance vs. exposure was measured for photochromic materials from 
the American Cyanamid Company, KBr crystals from Carson Laboratories and silver halide 
photochromic glass from the Corning Glass Works. When appropriate, the wavelength, tem- 
perature, and intensity were varied. In this section, typical results for KBr crystals and pho- 
tochromic glass are presented. The reciprocity effects will be discussed in section A. 7. 
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AMPLITUDE TRANSMITTANCE 


AKBr crystal is first colored to maximum density with UV radiation and then bleached 

o 

with 6328-A light derived from the He~Ne laser. Experiments were preformed in which a uni- 
form beam was allowed to strike a previously colored crystal. The intensity of the light passing 
through the crystal was measured with a photomultiplier tube and recorded as a function of time 
by means of a chart recorder. From these data, amplitude transmittance vs. exposure curves 
were obtained. Figure 27 shows the results for several crystals. 




EXPOSURE 


mW-sec 

2 

cm 


) 


EXPOSURE 


mW-sec 


cm 


(a) (b) 

FIGURE 27. TRANSMITTANCE vs. EXPOSURE CHARACTERISTICS OF KBr CRYSTALS. Constant temper- 
ature of 80°C. 


The effect of temperature on the sensitivity of the KBr crystals was determined by mea- 
suring amplitude transmittance vs, exposure while the crystal under test was held at a constant 
temperature with a thermoelectric heater unit. As shown in figure 28, the rate of bleaching 
increases as the temperature is raised, A temperature of 80°C is considered to be optimum 
by the manufacturers, since above 80°C, the thermal coloration effects tend to cancel any im- 
provement in bleaching sensitivity. 

Amplitude transmittance vs. exposure was also measured for samples of photochromic 
glass. Figure 29 shows results for three different wavelengths. All data were obtained at 
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FIGURE 28. TRANSMITTANCE vs. EXPOSURE CHARAC- 
TERISTICS OF KBr CRYSTALS SHOWING THE EFFECT 
OF TEMPERATURE VARIATIONS 


room temperature. The fact that this photochromic material is sensitive to a range of wave-* 
lengths indicates its application to multicolor holography (see sec. A.9). 

The exposure sensitivity data presented here show that for typical intensities used in 
2 

holography (<1 mW/cm ), the exposure times may be on the order of several minutes. The 
sensitivity of these materials is approximately 10 4 times less than Kodak 649F spectroscopic 
plates. 

The shape of the amplitude transmittance vs. exposure curves for KBr crystals can be 
justified analytically by assuming that the optical density D is directly proportional to the 
number of color centers N in the crystal,' That is, 

D = aN = log ^ 


75 




EXPOSURE 


mW-sec\ 

2 

cm j 


FIGURE 29. TRANSMITTANCE vs. EXPOSURE CHARACTERIS- 
TICS OF PHOTOCHROMIC GLASS SHOWING EFFECT OF WAVE- 
LENGTH VARIATIONS 


X » 

where =■ is the intensity transmission. The energy per unit area absorbed in the crystal 
in 

is ®in ~ ^out^ 111 turn ^ <l uantit Y proportional to the decrease in the number of color 
centers, 

(I. - I Jt = -rN 

' in out 

where r is a proportionality constant related to the bleaching efficiency. 

If we let I be the incident intensity, we have in the limit the following nonlinear differen- 
tial equation 

,f> . i (e - aN » . i) 

and the solution is 
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where is the initial number of color centers. 
Since 


T 1 

i alt 

1 - N ae r 
o 


the amplitude transmittance T a as a function of the exposure It is 



2 ^ 

The quantities N Q a and — can be considered constants depending on the particular crystal 
Figure 30 shows a typical curve which agrees quite closely with the experimental curves 
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cm 


FIGURE 30. CALCULATED TRANSMITTANCE vs. EXPOSURE 
CHARACTERISTICS OF KBr CRYSTALS 
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A. 3. ANGULAR ORIENTATION SENSITIVITY 

When the wavelength or the direction of incidence of the readout beam is varied contin- 
uously, there will be a corresponding change in the diffraction intensity from a hologram. The 
experimental methods used to determine the sensitivity of the hologram are either to vary the 
wavelength while the direction of incidence is fixed or to maintain a constant wavelength and 
vary the direction of incidence. The first method is difficult to implement because the wave- 
lengths of available sources cannot be varied continuously. The second and most convenient 
method can be carried out experimentally by rotating the incident beam relative to the holo- 
gram, or by rotating the hologram relative to the incident beam. It should be noted that the 
direction of incidence has two degrees of freedom. However, in all our experiments only one 
was employed. 

A.3,1. CALCULATED RESULTS. The angular width between half-power points of the 
curve of diffracted intensity vs. readout-beam angle can be derived from existing x-ray dif- 
fraction theory and is given by 

An _-\l & n 2 X 

where X is the constructing and reconstructing illumination wavelength 
6^ is the Bragg angle 

SD 

D is the effective thickness of the hologram 

If the incidence angles of the reference beam and signal beam are denoted by 9^ and 9 g 
(see fig. 31 ), each impinging from the two sides of the normal to the hologram, then 


Photo chromic Material 



Reference Beam 


FIGURE 31. DIAGRAM OF ANGLE 
CONVENTION FOR REFERENCE AND 
SIGNAL BEAMS 
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and 
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where t Q is the actual thickness of the hologram. 
The half -power width becomes 
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(60) 


This equation was obtained on the assumption that the absorption and the interaction between 
incident and diffracted waves in the hologram can be neglected. As we shall subsequently * 
demonstrate, the half -power width value in equation 60 is in good agreement with our experi- 
mental results. The approximate formula for is indeed applicable to a variety of angu- 

lar orientations and thicknesses of photo chromic materials. 

The equation for Afl.^ represents the internal half-power width. By applying Snell's re- 
lation, we can determine Afl^y^ in terms of the external parameters 
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which can also be represented as 
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In order to find the external half -power width value of the angular orientations (A^^air 7 
we reapply Snell's law to get 
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and finally 
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If the small-angie approximation of sin A 6 = A 9 is used, we then obtain the following 
relation 
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and in more complete form 
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This relation was used in all our calculations. Using equation 61, we have calculated the half- 
power width as a function of hologram thickness and as a function of reference beam and signal 
beam orientations. The results are shown in figures 32 through 35. The curves shown in fig- 
ure 28 depict the calculated external half -power width as a function of hologram thickness for 
three different offset angles. The offset angle, which as the angle subtended between the refer- 
ence beam and signal beam, was chosen so that its bisector is the normal to the hologram 
plane. The fringe surfaces were thus perpendicular to the hologram medium and spacing be- 
tween these fringe surfaces decreased with higher offset angle. In all instances, the Afl.^ * s 
inversely proportional to the thickness. 


The curve shown in figure 33 shows the half -power width as a function of incidence angle. 
Here 0^ = 0 g = 0 , so that the fringe surfaces are again perpendicular to the surface of the re- 
cording medium. 


The external half-power width as a function of the reference beam angle, for different sig- 
nal beam angles, is shown in figure 34. The thickness t Q represents a typical value and is a 
constant for these curves. Finally, we calculated the half -power width as a function of refer- 
ence beam angle, with the offset angle kept constant. The resulting data are shown in figure 35. 


As is evident from figures 33 and 34, there is an optimum value of 0^ for which the angu- 
lar orientation sensitivity is maximized. By differentiating ( A ^/ 2 ^air with res P ect anc * 

setting the result equal to zero, we have determined this optimum 0 when 0 = 0 ' and n and 

r r s 

0 g are fixed. For the first condition we obtain the relation 
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THICKNESS (fan) 

FIGURE 32. CALCULATED HALF-POWER WIDTH AS A 
FUNCTION OF HOLOGRAM THICKNESS FOR DIFFERENT 
OFFSET ANGLES, n =-1.49,X = 6328 A, 6 r = 


INCIDENCE ANGLE (deg) | 


FIGURE 33. CALCULATED HALF-POWER WIDTH 
AS A FUNCTION O^ INCIDENCE ANGLE, n = 1.49, 
X = 6328 A', e r = 6 S , t Q = 3175 /*. 



REFERENCE BEAM INCIDENCE ANGLE - 0 ; . (deg) 

FIGURE 34. CALCULATED HALF-POWER WIDTH 
AS A FUNCTION OF REFERENCE ANGLE FOR DIF- 
FERENT SIGNAL BEAM ANGLES, n = 1.49, X = 

6328 A, t 0 = 3175 fim. 



REFERENCE BEAM INCIDENCE 
ANGLE - S. (deg) 


FIGURE 35. CALCULATED HALF-POWER WIDTH 
AS A FUNCTION OF REFERENCE ANGLE WITH A 
FIXED OFFSET ANGLE OF 60°. t = 1587.5 #tm. 
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2 1 opt 

n = * 

2 sin 6 - 1 

r opt 

Using this relation, we calculated 8 as a function of the index of refraction n. A plot of the 

resulting data is shown in figure 36, The same results were also obtained when 0 - 0 o For 

s 

silver halide photochromic glass, n = 1.49 and consequently the optimum 0 is found to be 49.2°. 



INDEX OF REFRACTION - n 

FIGURE 36. CALCULATED OPTIMUM REFERENCE BEAM ANGLE AS A FUNCTION OF THE 
REFRACTIVE INDEX OF THE RECORDING MEDIUM. 0 r = f? s or 6> s = 0. 


For the more general condition of a fixed 8 , we obtain the following relation 
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Using this relation, we have determined 0 as a function of 6 for different fixed n. The 

r opt s 

resulting data were then plotted and are shown in figure 37* 

As is clearly indicated in our calculated results, the maximum angular orientation sen- 
sitivity, and consequently the least interaction between successive images, occurs at some 
optimum angle of incidence. This optimum angle is primarily a function of the index of refrac- 
tion of the hologram medium n and to some extent a function of signal beam angular orientation 



OPTIMUM 0 (degrees) 

FIGURE 37. CALCULATED OPTIMUM REFERENCE BEAM ANGLE 
AS A FUNCTION OF SIGNAL BEAM ANGLE, n = 1.49, n = 1.0. 

A. 3. 2. EXPERIMENTAL RESULTS. .All of our experiments involve the recording of sim- 
ple interference gratings which result from the interference of two plane waves. The effect of 
different hologram thicknesses and different offset angles between the reference beam and sig- 
nal beams was determined. The recording medium was either KBr crystals or silver halide 
glass. For comparison purposes, we also used the Kodak 649F spectroscopic plates which 
have a nominal emulsion thickness of 15 fim . 
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The experimental setup is shown in figure 38. The light from the laser is split into two 
paths by the beam splitter and subsequently diverted by the mirrors to the photochromic sam- 
ple. The pinhole assemblies and collimating lenses expand and collimate the beam. The re- 
sult is that two plane waves, I and I , interfere at the sample. Experiments were performed 

r s 

in which the interference pattern was recorded for various values of the angle 0 between the 
two beams and for various angular orientations of the photochromic material with respect to 
the light beams. Theorientation can be adjusted with great precision by means of a motor- 
driven platform. The angle 0 was changed by relocating the optical components. 



FIGURE 38. DIAGRAM OF SETUP FOR ANGULAR ORIENTATION SENSITIVITY 

MEASUREMENTS 

During readout, only the 1^ beam was allowed to be incident on the photochromic material. 
The intensity of the diffracted beam was measured with a photomultiplier tube. The resulting 
voltage was then fed into a paper chart recorder via a logarithmic amplifier. By rotating the 
photochromic at a constant rate, the relative intensity of the diffracted light was continuously 
plotted as a function of the incidence angle of the readout wave. 

It should be noted that the angle of diffraction 0 ^ must satisfy the grating equation 

sin 0 ^ + sin 6 ^ = ^ 
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where \ is the wavelength of the readout wave 
d is the fringe spacing of the grating 
0. is the angle of incidence of the readout wave 

i t 

0^ is the angle of the diffracted wave 

In our measurements, the readout source and the photomultiplier tube remained fixed 
while the three-dimensional grating was rotated. This means that the relation between & and 
8^ implied by the experimental configuration is that the sum of the incidence angle of the read 
out wave and the angle of diffraction is a constant. Figure 39 shows an example to illustrate 
this point. The curved line represents the actual variation based on the grating equation, 
whereas the straight line is the assumed variation. The difference between the two curves 
represents an experimental error. Since the curves coincide at the Bragg angle and remain 
extremely close for several degrees on either side, the error is negligible if measurements 
are confined to this region. Moreover, the orientation sensitivity of the thick recording mate 
rials under investigation here is so great that very little light is diffracted when the incidence 
angle of the readout wave differs from the Bragg angle by more than several minutes of arc. 



READOUT INCIDENCE ANGLE -8. (deg) 

FIGURE 39. CURVE REPRESENTING ERROR DUE TO FIXED PHOTOMULTIPLIER 

POSITION 
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Because of this, it was not necessary in these experiments to track the diffracted beam or use 
a collecting lens when measuring the intensity while the photo chromic sample was being ro- 
tated. 

Several examples of the experimental results are shown in figures 40, 41, and 42. Fig- 
ure 40 shows the relative diffracted intensity as a function of incidence angle of the readout 
beam for KBr crystals with two different thicknesses and offset angles. Similarly, figure 41 
shows results obtained from photochromic glass of three different thicknesses. Figures 41d, 
e, and f demonstrate the change in orientation sensitivity due to different reference beam 
angles even though the angle between the signal beam and reference beam is held constant. 

For comparison purposes, an orientation sensitivity curve obtained for Kodak 649F spectro- 
scopic plates is shown in figure 42. 

The half-power width of the angular orientation sensitivity curves was determined from 
the experimental data and compared with calculated results. These data are shown in table XI. 
Good agreement between calculated and experimental results is clearly evident. 

A.4. DIFFRACTION EFFICIENCY 

Diffraction efficiency is defined as the ratio of the diffracted light intensity to the incident 
light intensity. The experimental setup for measuring diffraction efficiency was essentially 
the same as that used in the orientation sensitivity experiments (see fig. 38). Simple inter - 

o 

ferograms of two monochromatic (632 8-A) plane waves were recorded in the photochromic 
materials under investigation. In the reconstruction process a plane wave was also used and 
the percentage of light diffracted by the photochromic was measured by means of a photomul- 
tiplier tube. 

Efficiencies were measured for interference gratings made with different exposure times, 

reference-beam-to-signal-beam ratios, and different spatial frequencies of the fringes as 

determined by the angle between the two beams. 'Typical results obtained for KBr crystals 

range from 0.002% to 0.175% for exposures of 480 mW ~^ ec and 2400 - ~^ ec , respectively, 

crn cm 

with an offset angle of 40°. Figure 43 shows the variation in diffraction efficiency as a func- 
tion of exposure with relatively low exposures of photochromic glass. 

If an interferogram is recorded in a photochromic material and then reilluminated with 
one of the beams of light, the diffraction efficiency will decrease as shown in figure 44. This 
behavior has important consequences in the multiple -image storage problem which is dis- 
cussed in section A. 5. 

A.5. MULTIPLE-IMAGE STORAGE 

One of the more important reasons for using a three-dimensional recording medium is 
its potential application to the storage of multiple signals. In this section we will discuss the 
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CHANGE IN INCIDENT ANGLE OF READOUT BEAM (min) 

(a) Thickness = 4.25 mm. Offset angle = 60°, 8 r - (9 g = 30°. 


m 



CHANGE IN INCIDENT ANGLE OF READOUT BEAM (min) 

(b) Thickness = 2 mm. Offset angle = 60°, 0 r = 0 g = 30°. 



CHANGE IN INCIDENT ANGLE OF READOUT BEAM (min) 

f 

(c) Thickness = 2 mm. Offset angle - 40°, 6 = 0 = 20°. 

r s 


FIGURE 40. EXPERIMENTAL RESULTS OF ANGULAR ORIENTATION 
SENSITIVITY MEASUREMENTS OF KBr CRYSTALS 


8 ? 


CD 



(a) Thickness = 1587.5 ix m. Offset angle = 10°, 6 r = G s ~ 5 °* 



(b) Thickness = 3175 ^m. Offset angle = 10°, 6^ - 0 g = 5°. 



(c) Thickness - 6350 jum. Offset angle - 10°, = 5°. 


FIGURE 41. EXPERIMENTAL RESULTS OF ANGULAR ORIENTATION SENSITIVITY MEASUREMENTS OF 

SILVER HALIDE GLASS 
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CHANGE IN INCIDENT ANGLE OF READOUT WAVE (min) 

(d) Thickness = 1587.5 Offset angle = 60°, 0 = 0°, $ = 60°. 



CHANGE IN INCIDENT ANGLE OF READOUT 
WAVE (min) 


(e) Thickness = 1587.5 fim. Offset angle = 60°, Q = 30°, 0 = 30°. 

r s 



CHANGE IN INCIDENT ANGLE OF READOUT WAVE (min) 


(f) Thickness = 1587.5 juni. Offset angle = 60°, 0 = 60°, 0 = 0°. 

» r s 

FIGURE 41 EXPERIMENTAL RESULTS OF ANGULAR ORIENTATION SENSITIVITY MEASUREMENTS OF 

SILVER HALIDE GLASS (Continued) 
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CHANGE IN INCIDENT ANGLE OF READOUT 
WAVE (min) 

FIGURE 42. EXPERIMENTAL RESULTS OF ANGULAR 
ORIENTATION.SENSITIVTTY MEASUREMENTS OF KODAK 
649-F SPECTROSCOPIC PLATES. Offset angle = 60°, $ v = 

e s = so 0 . 


TABLE n. COMPARISON OF CALCULATED VALUES OF HALF-POWER WIDTH 
OF ANGULAR ORIENTATION SENSITIVITY CURVE WITH EXPERIMENTAL VALUES 


Reference Beam 
Angle 9 r 

(deg) 

Signal Beam 
Angle 

(deg) 

Thickness 

(fim) 

Calculated Half- 
Power Width 
(min) 

Experimental 
Value of Half- 
Power Width 
(min) 

60 

0 

1587.5 

3.6 

3.2 

30 

30 

1587.5 

2.1 

2.0 

0 

60 

1587.5 

1.8 

1.7 

5 

5 

1587.5 

11.0 

10.9 

5 

5 

3175.0 

5.5 

5.3 

5 

5 

6350.0 

2.75 

2.8 
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EXPOSURE 



EXPOSURE 



FIGURE 43. CURVE OF DIFFRACTION EFFICIENCY AS 
A FUNCTION OF EXPOSURE LEVELS FOR SILVER HALIDE 
GLASS. Thickness = 3175 /z, I r = I g = 0.2 mW/cm 2 , 8 t = 

e s = 5°. 


FIGURE 44. CURVE ILLUSTRATING REDUCTION OF 
DIFFRACTION EFFICIENCY DURING RECONSTRUC- 
TION. Thickness = 3175 P-, I r = I s = 0.27 mW/cm 2 , 

8 r = $ s = 5°, exposure time = 10 min. 


practical aspects of the photochromic devices in actual system operation and the problems 
encountered. 

The Bragg discrimination must be considered in relation to limitations imposed on the 
number of holographically stored signals, when we store many complex images each con- 
structed with a differently oriented reference beam. With photographic plates having an emul- 
sion thickness of 15 jam, the average change of orientation between each successive image with 
little interaction is 10°. The angular orientation sensitivity is quite low at low offset angles 
and increases with higher offset angles. Since the available angular orientation in constructing 
a conventional hologram, is 180°, we are confined to storage of only 18 images. With special 
coding techniques, separation of images, and wavelength variations, the calculated number of 
stored images can be significantly increased. In this section we deal solely with Bragg 
holography. 

A calculated increase of two to three orders of magnitude can be achieved in stored images 
with existing photochromic materials. In practice, however, the primary limitation is not the 
number of images. The diffracted intensities of recordings made early in a sequence of equal 
recordings are either below the noise level or at best very weak. This anomalous behavior, 
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which is peculiar to photochromic materials and is not evident in photographic material, is 
believed to be caused by some nonlinear bleaching effect. Assuming that the photochromic 
material is comprised of a large number of layers of equal thickness, the layer closest to the 
illumination source is partially bleached with the first recorded image. Illumination of subse- 
quent images continues to bleach the first layer and penetrates into the inner layers. The 
cumulative bleaching of the first layer continues until the saturation point when the recorded 
information about the first image is obliterated. This process is then repeated at the second 
layer, at the third layer, and so on, each time obliterating the recorded information or re- 
ducing the recorded fringe contrast and thereby lowering the diffraction. efficiency. We have 
performed a number of experiments which demonstrate the nonlinear bleaching behavior in 
photochromic glass and KBr crystals and have found that the latter create the most difficult 
problems. 

« Early experiments involved the recording of simple interference patterns of two plane 
waves, each at a different orientation of the photochromic material. To date a sequence of 40 
separate gratings, each incoherent with the others, has been successfully recorded in the KBr 
crystals, whereas around 200 gratings were recorded in the photochromic glasses. Two ex- 
amples of these results which illustrate the cumulative bleaching problem are shown in fig- 
ure 45 and figure 46. The curves represent fluctuations in diffracted intensity as the orienta- 



FIGURE 45. EFFECT OF NONLINEAR BLEACHING ON KBr CRYSTALS 



FIGURE 46. EFFECT OF NONLINEAR BLEACHING ON SILVER HALIDE GLASS 
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tlon o f the incident illumination with respect to the photochromic medium is continuously altered. 
Tlie angular separation between recordings was sufficiently large so that the SJjLi structure 
for some recordings is clearly evident. 

For the examples shown In figures 45 and 46, 10 recordings were made In each medium* 

The average offset angle was approximately 10° and the change in angular orientation lie tween 

each recording was approximately four minutes of arc* The intensities of the reference beam 
2 

and signal beam were 2 mW/cm and the exposure time for each recording was five minutes* 

For the KBr crystal, the effect of nonlinear bleaching is dramatically demonstrated in fig- 
ure 45. The diffracted intensity of the tenth recording is shown at the extreme left and Is the 
highest* The diffracted intensities of the earlier recordings decrease down to the sixth re- 
cording, and the remainder cannot be detected. The situation Is improved for the photochromic 
glass, and in figure 46 the diffracted intensity variations for all lt> recordings can be seen 
clearly. The diffracted intensity of the tenth recording is on the extreme right and although 
there is some effect of nonlinear bleaching, it is not significant* 

We also recorded 100 relatively complex diffraction patterns in one hologram on the l/g- 
inch -thick photochromic glass. The jatterns resulted from the Interference between plane 
waves of the reference beam and complex waves from the signal beam, A sequence of 100 
diffusely muminated transparencies were used as objects, each having a bandwidth of approx - 
imately 150 lines, mm. Exposure times for each recording ranged from two to three minutes, 
and angular orientation (between reference beam and hologram) was changed by four minutes 
of arc to separate each recording. 

Finally, the reconstructed image of each recording was investigated. A representative 
sample of the results Is shown in figure 47. As expected, the signal -to -noise ratio for the 
later recordings was relatively high, providing reconstructed images of high quality* As is 
evident from the results shown in figure 47a and b, the signal -to -noise ratio of the recon- 
structed images is significantly degraded in the earlier recordings, but in no instance were 
the images totally obliterated. It is believed that by optimizing the illumination intensities and 
exposure time, 200 or 300 complex recordings can easily and effectively be stored in existing 
materials* 

An illustration of the excellent Bragg discrimination capabilities for photochromic glass 
is shown in figure 48. There Is very little if any, interaction between recordings. The ob- 
jects chosen were transparencies showing the face of a time clock, The first recording was 
made at 0 minutes and 0 seconds and In subsequent recordings the second hand was advanced 
by 15 seconds* On readout, if the hologram Is given continuous angular motion, the recon- 
structed image of the clock face appears to have Its second hand move continuously. This, in 
effect, is a stored movie* The 100 recordings In the photochromic material replace 100 
frames of strip film. 
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(a) Imago from Fifth Recording 


(bj Image from Twentieth Recording 



FIGURE 4?. RECONSTRUCTED IMAGES FROM A MULTIPLE- EXPOSURE HOLOGRAM OF 
100 RECORDINGS. Angular separation between each recorded Image la 4 min of arc. 
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fel Imago from Ninetieth Recording 

FIGURE 41 RECONSTRUCTED IMAGES FROM A MULTIPLE -EXPOSURE HOLOGRAM OF 
100 RECORDINGS (Continued > 


A.6, SHRINKAGE EFFECTS 

The presence of shrinkage In three-dimensional recording media generally causes a 
change in both orientation and spacing oi the fringe surfaces in the hologram. As a consequence 
of this change, either the orientation or the wavelength of the reconstructing beam must be al- 
tered to maximize the diffracted intensity. In both instances, the location of the reconstructed 
image will also change. In photographic material, the shrinkage is caused by the emulsion 
shrinkage that occurs during chemical processing. In photochromlc crystal, the shrinkage is 
due to reduction of temperature between readin and readout. 

The shrinkage phenomenon inKBr crystals was investigated by calculating the effect of 
temperature reduction in the crystals on the orientation of the reconstructing beam and con- 
firmed our results experimentally. A simple grating hologram formed by the interference of 
two coherent plane waves was used. In three-dimensional media, the recorded interference 
pattern is composed of parallel plane surfaces. The orientation and frequency of these sur- 
faces are described in relation to the schematic illustration of figure 49. For simplification, 
symmetrical orientation for the two plane waves was chosen which produced fringe surfaces 
perpendicular to the surface of the recording medium. The spacing of the recorded fringes d 
obeys the Bragg relation 

d -FlfiTeT 
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(c) imags from Sixty-third He carding (d) Image from Sixty-fourth He cording 


FIGURE 48. RECONSTRUCTED IMAGES FROM A MULTIPLE -EXPOSURE HOLOGRAM OF 
100 RECORDINGS ILLUSTRATING BRAGG DISCRIMINATION CAPABILITIES Angular sepa- 
ration between each recorded image la 4 min el arc. 
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[el Image from Sixty-fifth Record Lng 

FIGURE 48. RECONSTRUCTED IMAGES FROM A MULTIPLE -EXPOSURE HOLOGRAM OF 
100 RECORDINGS ILLUSTRATING BRAGG DISCHIMIN ATI0 N CAPABILITIES (Continued) 




FIGURE 4&. DIAGRAM ILLUSTRATING FRINGE SURFACES 
AND ORIENTATION OF WAVES 
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where X is the illumination wavelength, and tMs hail the angle subtended between the two plane 
waves. 

Using the thermal expansion coefficient of 43 * 10~®/°C lor KBr and denoting the temper- 
ature difference between read in and readout as AT, we obtain an expression for the fringe 
spacing alter contraction 

The angle satisfying the Bragg relation lor the contracted case is 

. a , sin Q 

am w — — t’ — 

1 - 43 x io AT 

and the required angular reorientation between rend in and readout becomes 

&* - 6 = sin’ 1 ( sin# \ _ 6 f63) 

\l - 43 x iff* AT/ 

It is interesting to note that the required, reorientation is not a function of the thickness of the 
crystal. 

In our experiment we used the values 


and 


a = 

Fi _ = so^c 


T , . = £5°C 

readout 


AT 


c 


^ 55 °C 


By substituting these parameters into equation 53, we find that 
9* - 9 - 4,6 min of arc 
which is in exact agreement with the experimental result. 

Although the required reorientation is not a function of crystal thickness, the thickness 
parameter plays an important role. This Is particularly evident when complex gratings of a 
variety of fringe spacings are recorded. In such situations, the shrinkage effects on fringe 
spacing are varied so that the Bragg relation cannot be simultaneously satisfied lor all fringes. 
The deleterious effect of nonuniform image Intensity can then be observed. As the thickness 
of the crystal is increased, this problem Is more severe because of greater angular orientation 
sensitivity. The photographs shown In figure 50 show the undesirable effects of a reduction in 




fa) Temperature 90°C 


(c) Temperature 70°C 


fb) Temperature 75°C 


ft)) Temperature 65°C 


FIGURE 50. DEGRADATION OF RECONSTRUCTED IMAGE DUE TO REDUCTION IN TEMPERATURE 




temperature of a crystal during the reconstruction process. If at any temperature the holo- 
gram is reoriented by the proper amount, the quality of the reconstruction is essentially the 
same as the reconstruction obtained at 8fl°C, 

A new technique proposed by Carson Laboratories can eliminate the shrinkage problem in 
KBr crystals. Retardation by red light or coloration of crystal is used. The recording takes 
place at room temperature using the He-Ne laser, while the crystal is being colored with in- 
coherent UV radiation. This simultaneous illumination by UV and red light Is essentially con- 
trolled coloration Of the crystal* Shrinkage is eliminated, since it is no longer necessary to 
record at elevated temperatures and reconstruct at lower temperatures* In the phoEochromic 
glass and in some of the photochromic plastics, the problem of shrinkage is nonexistent since 
both readin and readout are done at room temperature* 

A,7, RECIPROCITY FAILURE 

Reciprocity failure provides a measurable criterion for determining the inefficiency of a 
recording medium at different levels of intensity relative to Its efficiency at the optimum in- 
tensity* This measure plays an important role in wavefront reconstruction, particularly when 
a wide range ol Intensity levels are to be recorded* 

A series of experiments was performed to measure the reciprocity failure of the KBr 
crystals* The amplitude transmittance was measured as a function of exposure for a 50- dB 
range of incident light intensities. The results for two different crystals are shown in Fig- 
ure 51, The reciprocity curves In each figure were measured for different bias levels of am- 
plitude transmittance* In aU cases, the results show severe high -intensity reciprocity failure. 
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Maximum Optical Density ■ 1*7 


FIGURE 51* RECIPROCITY CURVES FOR KBr CRYSTALS. Thickness - 2 mm. 
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The implication of these results Is indeed discouraging. We had hoped that by resorting to 
very intense illumination sources such as pulse lasers, we could appreciably reduce the readin 
time. It is evident from the results that this did not happen. The reciprocity law does rtot 
hold at any ol the Investigated intensity levels, 

A .8, RECORDING COMPLEX INTERFERENCE PATTERNS 

In order to demonstrate the use of photochromlc materials in holography, experiments 
were performed in which the signal beam was spatially modulated by transparencies using both 
planar and diffuse illumination. Two different experimental setups for recording were used. 

In one, both the reference beam and signal beam were introduced from the same side of the 
recording medium. In the other, the reference beam and signal beam were introduced from 
opposite sides. The readout orientation sensitivity is relatively broad for the interference 
pattern resulting from the latter setup. The spectral resolution, on the other hand, is extreme- 
ly narrow since the recorded gratings act as spectral interference filters. Readout with a 
mercury -arc source or a zirconium lamp produced single -color diffracted images of narrow 
spectral bandwidth. 

It is interesting to note that the spectral sensitivity of thick photochromlc materials is suf- 
ficiently great even for low offset angles that a back-reflection type of hologram is not required 
in order to obtain good reconstructions with white light. Figure 52 shows the reconstructed 
image Of a signal comprised of transparent targets on a dark background. For the recording, 
the signal was coherently illuminated, and both the reference beam and the signal beam were 
introduced from the same side of the crystal. Figure 53 shows the reconstructed image of a 



FTGUR E 52 . RECONSTft UCT ED IMA C E OF COM ERENTLY ILLU- 
MINATED TRANSPARENT OBJECTS ON A DARK BACKGROUND 
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FIGURE 53, RECONSTRUCTED IMAGE OF A DIFFUSELY ILLU- 
MINATED, TRANSPARENT* CONTINUOUS-TONE OBJECT 

continuous -tone transparent signal which was diffusely illuminated during recording. In this 
case, the reference beam and the signal beam were introduced from opposite sides of the re- 
cording medium* The recording and readout were clone with the 6 32 8- A illumination from a 
He - Me laser* In both cases, the quality ol the reconstructed images was comparable to re- 
suits with Kodak 649 F spectroscopic plates. 

A . 9. MU LTlC OLQR HGLOG KAPH Y 

As shown earlier, the silver halide photochromic glass is sensitive to various illumination 
wavelengths and therefore can be effectively employed in multicolor wavefront reconstruction. 
A number of experiments using photochromic glass was performed* A hologram made of 
three superimposed diffraction patterns was constructed. Each diffraction pattern resulted 
from the interference of reference and signal waves of only one of three Illumination wave- 
lengths, which were the 6328-A wavelength derived from a He-Ne laser, and the 4880-A and 
5145- A wavelengths from the argon laser. The signal used in our experiments was a diffusely 
illuminated multicolor transparency and the reference beam was a plane wave. The average 
angle subtended between the reference beam and signal beam was approximately UF. The 
construction geometry was such that the reference beam for all three wavelengths was derived 
from a common point source* 

The three wavelengths were then used in the reconstruction, again derived from a common 
point source, and the resulting reconstructed image was investigated* The reconstructed 
image, shown in figure 54 is ol high quality and duplicates the color hues of the actual signal. 

In addition, the ghost images and attendant detrimental effects were completely absent. For 
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comparison, we used the same experimental setup and constructed a hologram with Kodak 649 F 
spectroscopic plates. As shown In figure 55, the reconstructed image contains six ghosE 
images. Some of them overlap on the desired image, resulting in a relatively poor reconstruc- 
tion. 

As demonstrated by the experimental results, there are obvious advantages in using the 
thick photochromtc glass instead ol photographic film* Unfortunately, the exposure times in- 
volved In our experiments were IQ 4 times longer for the photochromic material than lor the 
film. 



FIGURE 54 RECONSTRUCTED IMAGE OF A DIF- 
FUSELY ILLUMINATED MULTICOLOR TRANS- 
PARENCY USING SILVER HAUDE GLASS. Thickness 
= 3175 Mm, offset angle * 16°. 



FIGURE $5. RECONSTRUCTED IMAGE OF A DIF- 
FUSELY ILLUMINATED MULTICOLOR TRANS- 
PARENCY USING KODAK 4H9-F SPECTROSCOPIC 
PLATES. Thickness ■ Is pm, offset angle = 10°. 
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